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1 .  Introduction 


Direct  measurements  of  the  infrared  radiance  from  the  upper  atmosphere  and, 
especially,  of  the  variability  in  this  radiance  are  very  limited.  Furthermore, 
almost  no  data  exists  on  the  variations  of  the  infrared  emissions  with  geophysical 
parameters,  such  as  latitude  or  magnetic  activity,  or  with  localized  dynamic 
effects.  Thus,  to  obtain  estimates  of  the  infrared  radiance  of  the  atmosphere, 
one  must  rely  on  computer  models.  The  computer  model  currently  being  used  to 
predict  the  infrared  radiance  at  high  altitude  is  a  modified  version  of  the  one  de¬ 
scribed  by  Degges  and  Smith^  (1977).  This  model  is  intended  to  simulate  the 
radiative  properties  for  average  midlatitude  atmospheric  conditions.  In  this  re¬ 
port  we  discuss  data  and  supporting  phenomenology  that  are  required  to  upgrade 
this  model  to  include  the  variations  of  the  upper  atmospheric  radiation  resulting 
from  changes  in  geophysical  conditions. 

The  intensity  of  infrared  emissions  in  the  upper  atmosphere  is  dependent  on 
both  the  density  and  distribution  of  specific  molecules  and  the  sources  that  excite 
the  emitting  levels  of  these  molecules.  At  the  present  time  there  is  insufficient 
data  concerning  the  density  and  the  distribution  of  the  IH  emitting  species  of  the 
upper  atmosphere  to  be  able  to  fully  describe  these  parameters.  However,  the 

(Received  for  publication  22  May  1981) 

1.  Degges,  T.  C.  and  Smith,  H.J.P.  (1977)  A  High  Infrared  Radiance  Model, 
AFGL-TR-77-0271,  AD  A059  242.  Final  Report,  Contract  F19628-77-C- 
0041,  Visidyne,  Inc.,  Burlington,  Maes. 
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density  and  distribution  of  these  species  can  be  estimated  with  varying  degrees  of 
confidence.  These  estimates  are  based  on  calculations,  which  are  supported  by 
measurements,  and  on  models  that  include  the  current  understanding  of  the  photo¬ 
chemical,  excitation,  and  transport  processes  in  the  upper  atmosphere.  In  this 
report,  data  that  can  be  related  to  either  the  infrared  emission  or  to  the  species 
concentrations  that  produce  these  emissions  are  considered. 

The  long-term  purpose  of  this  study  is  to  develop  a  data  base  sufficient  for 
the  modeling,  not  only  of  the  IR  radiance,  but  of  its  variability.  Both  spatial  and 
temporal  inhomogeneities  and  their  magnitudes  and  frequencies  of  occurrence  are 
considered.  This  report  represents  an  effort  to  clarify  and  summarize  the  state 
of  our  knowledge  at  the  present  time.  An  attempt  has  been  made  to  identify  those 
aspects  of  the  problem  where  additional  effort  needs  to  be  made. 

The  variations  and  structured  features  observed  in  the  concentrations  of  the 
upper  atmospheric  species  may  be  divided  into  three  general  classes.  The  first 
general  class  is  the  long-term  variations  of  the  quiescent  atmosphere  that  are 
associated  with  vertical  profiles.  These  variations  can  be  described  by  the  hydro¬ 
static  equation  and  the  ideal  gas  law.  From  a  knowledge  of  the  density,  temper¬ 
ature,  and  composition  near  the  earth's  surface,  the  profiles  of  the  major  species 
can  be  derived  from  measured  profiles  of  temperature,  turbulent  and  molecular 
diffusion  velocities,  solar  radiation  dissociation  rates,  and  chemical  reaction 
rates  (for  example,  see  Keneshea  et  al,  1979  ).  This  class  of  atmospheric  vari¬ 
ations  includes  those  resulting  from  atmospheric  heating  (which  depends,  in  turn, 
upon  solar  cycle,  latitude,  season,  and  local  time).  The  profiles  of  this  class  of 
variations  are  described  in  detail  by  well-developed  empirical  models  such  as 
CIRA  72,  USSA  76,  Jacchia  77  (J77).  and  MSIS  78. 

The  second  and  third  classes  of  atmospheric  variations  are  shorter  period 
waves  and  localized  heating  events  that  are  superimposed  on  the  smooth  diffusive 
profiles  of  the  first  class.  The  second  class  includes  those  short-term  wave-like 
structures  that  are  the  fluid  responses  to  the  many  different  types  of  forces  applied 
to  the  atmosphere.  These  result  in  structure  with  a  wide  range  of  scale  sizes, 
amplitudes,  and  irregular  features  that  result  from  wave  interference  or  dissipa¬ 
tive  breakdown  of  waves  into  turbulent  eddies.  The  source  of  these  variations  in¬ 
clude  the  diurnal  and  semidiui’nal  (and  higher  order)  tidal  motion  (Chapman  and 
Lindzen.  19711,  Forbes  and  Garrett,  1978  )  and  the  gravity  wave  structures  that 

2.  Keneshea,  T.J.,  Zimmerman,  S.P.  and  Philbrick,  C.R.  (197  9)  A  dynamic 

model  of  the  mesosphere  and  lower  thermosphere,  Planet.  Space  Sci. 

22:385-401. 

3.  Chapman,  S.  and  Lindzen,  R.S.  (1970)  Atmospheric  Tides,  D.  Reidel, 

Hingham,  Mass. 

4.  Forbes,  J.M.  and  Garrett,  H.B.  (1978)  Solar  diurnal  tide  in  the  thermosphere, 

J.  Atmos.  Sci.  35:148-208. 


are  generally  propagated  from  the  lower  atmosphere  into  the  middle  and  upper 
atmosphere  (Hines,  1974,^  Reber  et  al,  1975®). 

The  third  class  of  variations  includes  irregular  structures  that  are  due  to  the 
deposition  of  coupled  energy  from  magnetospheric  particles  and  fields.  These 
energy  sources  provide  strong  localized  heating  and  thereby  produce  many  types 
of  irregular  features  in  the  atmosphere.  Such  heating  may  result  from  direct 
particle  precipitation  into  the  atmosphere  from  the  solar  wind  or  the  trapped  par¬ 
ticle  regions  of  the  magnetosphere,  such  as  aurora.  However,  the  largest  heat¬ 
ing  effects  result  from  plasma  interactions  with  the  neutral  atmosphere  when 

n 

electric  and  magnetic  fields  accelerate  the  ionospheric  plasma  (Cole,  197  1,  Mayr 

O 

and  Harris,  1978  ).  The  high-latitude  regions  aie  frequently  subjected  to  large- 
scale  Joule  heating  and  to  very  localized  effects,  such  as  the  "momentum  source" 
and  high  structured  particle  precipitation  events.  The  strongly  heated  regions 
may  be  so  perturbed  by  the  energy  input  that  the  vertical  profiles  of  individual 
species  can  no  longer  be  represented  by  hydrostatic  profiles. 

The  modeling  efforts  and  experimental  data  will  be  reviewed  in  order,  accord¬ 
ing  to  the  three  classes  of  variations  described  above.  This  ordering  has  the  ad¬ 
vantage  of  allowing  a  continuity  of  development  in  the  complexity  of  the  atmospheric 
features.  In  addition,  the  variations  resulting  from  each  class  are  generally  con¬ 
sidered  a  superposition  to  the  effect  of  the  previous  classes.  Recently  a  great 
deal  of  attention  has  been  given  to  determining  the  atmospheric  properties  at  higher 
satellite  altitudes.  However,  the  regions  that  contribute  most  significantly  to  the 
infrared  emission  of  interest  are  from  50  to  200  km.  This  review  will  concentrate 
on  satellite  and  rocket  results  at  these  altitudes. 

The  application  of  the  atmospheric  models  and  measurements  to  the  calculation 
of  infrared  radiance  will  generally  take  two  forms.  In  the  first  form,  that  applies 
generally  to  the  smooth  variations  of  the  first  class,  an  atmospheric  profile  of 
temperature  and  species  concentration  is  used  as  an  input  to  a  model  for  the  in¬ 
frared  radiance  (such  as  the  High  Altitude  Infrared  Radiance  Model  of  Degges^ 
(Degges  and  Smith,  1977)  described  in  Chapter  2)  and  the  infrared  radiance  evalu¬ 
ated.  Examples  of  infrared  radiances,  based  on  calculations  using  the  atmospheric 


5.  Hines,  C.O.  (1974)  The  Upper  Atmosphere  in  Motion,  American  Geophysical 

Union  Monograph  Ser  18. 

6.  Reber,  C.A.,  Hcdin,  A.E.,Pelz,  D.  T.  Potter,  W.  E.  and  Brace,  L.  H.  (1975) 

Phase  and  amplitude  relationships  of  wave  structure  observed  in  the  lower 
thermosphere,  J.  Geophys,  Res.  80:4.576-4580. 

7.  Cole,  K.D.  (1971)  Electrodynamic  heating  and  movement  of  the  thermosphere. 

Planet  Space  Sci.  19:59-75. 

8.  Mayr,  H.G.  and  Harris,  I.  (1978)  Some  characteristics  of  electric  field 

momentum  coupling  with  the  neutral  atmosphere,  J.  Geophys.  Res. 
83:3327-3336. 
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properties  described  in  Chapter  3.  are  given  in  Chapter  4.  In  Chapter  3  var  ious 
atmospheric  models  are  discussed  with  a  concentration  on  the  J77  model.  Ex¬ 
amples  of  the  predictions  of  the  J77  model  are  presented  and  compared  to  meas¬ 
urements  of  the  atmospheric  properties.  In  addition,  existing  data  concerning 
the  variability  of  the  minor  atmospheric  species  is  also  reviewed. 

In  the  second  form,  where  irregular  variations  are  considered,  the  standard 
deviation  and  the  auto -correlation  function  or  power  spectral  density  are  used  to 
estimate  the  expected  variations  in  the  radiance  as  discussed  in  Chapter  5.  These 
estimates  require  certain  assumptions  concerning  the  connection  between  atmos¬ 
pheric  properties  and  the  resultant  radiance.  In  Chapter  5  the  sources  of  the 
small-scale  variability  are  discussed,  experimental  data  reviewed,  and  a  prelim¬ 
inary  summary  of  the  sizes  and  magnitude  of  the  variability  presented.  In  addition, 
an  example  of  the  use  of  this  material  to  construct  a  power  spectral  density  is 
given,  as  well  as  a  proposed  method  of  generating  a  model  for  the  auroral  power 
spectral  density  for  two  specified  wavelengths  in  the  infrared. 

Finally,  in  Chapter  6,  a  brief  summary  of  the  important  problems  that  must 
be  given  further  attention  and  that  point  the  directions  for  future  research  are 
briefly  discussed. 
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2.  Current  High  Altitude 
Infared  Radiance  Model 
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The  High  Altitude  Infrared  Radiance  Model  represents  the  ouri'cnt  capability 

1  3  ' 

in  a  conMnuing  program  (Corbin  tt  al,  19G0,  Dcgges,  1972,"  1974,  Degges  and 
Smith,  1977  )  to  predict  infrared  radiance  levels  in  the  upper  atmosphere  for  the 
spectral  range  from  2.5  to  2  5  pm.  In  these  predictions  emphasis  is  placed  on 
radiation  originating  at  altitudes  between  GO  and  500  km.  Comparison  of  the  pre¬ 
dictions  of  the  computer  model  with  experimental  measurements  serves  to  chei  k 
the  adequacy  of  the  infrared  model  and  provides  greater  understanding  of  the 
atmospheric  properties  that  determine  the  infrared  radiation. 


1.  Corbin,  V .  L. ,  Dalgarno,  A.,  Degges,  T.  C. ,  House,  F.  B.  ,  Lilienfeld,  P.  , 

Ohring,  G. ,  and  Oppel,  G.E.  (1969)  Atnmspheric  Radiance  Models  for 
Limb-Viewing  Geometry  in  the  Flve-to-Twenty-five  Micron  Spectral 
Region,  AFCRL-G9-0552,  AD  868  194L. 

2.  Degges,  T.  C.  (1972)  A  High  Altitude  Radiance  Model,  AFCRL-72 -0273, 

AD  745  319,  Final  Report,  Contract  F19628-71-C-0156,  Visidyne,  Inc., 
Burlington,  Mass. 

3.  Degges,  T.C.  (1974)  A  High  Altitude  Radiance  Model.  AFCRL-TR-74-0606. 

ADA008  035,  Final  Report,  Visidyne,  Inc.,  Burlington,  Mass. 

4.  Degges,  T.C.  and  Smith,  H.  J.  P.  (1977)  A  High  Infrared  Radiance  Model, 

AFGL-TR-77-0271,  Final  Report,  Contract  19628-77 -C-064l,  Visidyne, 
Inc.,  Burlington,  Mass.,  AD  A059242. 


The  principal  sources  of  infrared  radiation  at  altitudes  greatt  r  than  rm  kjn 
arc  water  vapor  (11., O),  carbon  dioxide  (CO.,),  ozone  (Hy).  nitrii-  o\ult  (NO),  .  ar- 
bon  monoxide  (CO),  and  the  hydroxyl  i-adical  (OH).  The  i-esults  of  the  recent 
SPIKE  experiment  (Stair  et  al,  i;i81'’)  are  consistent  with  this  conclusion.  W  hile 
some  calculations  have  been  made  that  include  CO;  its  abundance  and  altitude 
variation  are  uncertain.  It  is  not  included  in  the  model.  Similaily,  while  OH  is 
known  to  be  important  at  wavelengths  from  the  visible  through  3.  3  nni,  it  is  not 
presently  included.  Although  the  model  includes  altitude  profiles  of  the  species 
concentrations  of  H.,0,  CO.,,  O^,  and  NO,  there  is  not  suffic  ient  knowledge  to 
provide  con.  ontration  profiles  for  all  latitudes  and  seasons.  The  present  model 
approximates  average  midlatitude  conditions. 

The  radiance  model  is  further  limitt'd  by  the  incomplete  knowledge  of  the 
physical  processes  that  control  the  population  of  infrared  emitting  states  of 
atmospheric  molecules.  Except  for  the  pure  rotational  radiation  from  molecules 
such  as  water,  the  degree  of  excitation  of  vibrational  le  vels  determines  the  radi¬ 
ation  from  infrared  emitting  molecules.  The  most  important  mechanisms  for  ex¬ 
citation  of  vibrational  and  rotational  levels  are  collisional  excitation  and  deexcitation, 
absorption  and  reemission  of  electromagnetic  radiation,  and  exothermic  chemi- 
lumiescent  reactions.  In  the  troposphere  and  lower  stratosphere  collisional  pro¬ 
cesses  are  rapid  enough  to  control  the  population  of  vibrational  levels.  Above 
30  km,  however,  collisional  excitation  bec;omes  less  efficient  and  radiative  pro¬ 
cesses  become  important  so  that  the  combined  effects  of  collisional  and  radiative 
processes  must  be  considered.  Below  an  altitude  of  about  PO  km,  molecular 
nitrogen  and  oxygen  arc  the  most  important  collision  partners.  .Xbove  this  altitude 
atomic  oxygen  becomes  important  in  determining  the  collisional  excitation  and 
deexcitation  rates. 

In  determining  the  effects  of  radiation  on  the  populations  of  the  vibrational 
levels  of  infrared  active  molecules,  it  is  necessary  to  separate  the  radiation  due 
to  a  single  change  in  vibration  quantum  numbers  from  the  rest  of  the  radiation 
field.  This  can  be  achieved  by  assuming  that  the  individual  rotational  lines  of  a 
band  have  Doppler  spectral  profiles  and  do  not  overlap.  This  assumption  allows 
the  construction  of  radiative  transfer  functions  that  are  adequate  numerical 
approximations  for  the  calculation  of  the  absorption  of  infrared  energy  from  the 
sun,  the  earth's  surface,  and  molecules  at  all  altitudes  in  the  atmosphere. 

Figure  1  illustrates  the  geometry  employed  in  the  radiance  model.  A  molecule 
at  point  P  absorbs  radiation  from  three  sources:  1)  the  sun,  2)  the  lower  boundary 


5.  Stair,  A.  T.  ,  Jr.,  Nadile,  R.M.,  Ulwick,  J.  C.  ,  Baker,  K.  D.  and  Baker,  D.J. 
(1981)  Infrared  measurements  of  aurora  airglow  and  the  upper  atmosphere. 
Paper  presented  at  AIAA  Aerospace  Sciences  Meeting,  Jan.  1981,  St.  Louis, 
Missouri,  Paper  no.  AIAA-81-0421. 
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F  igure  1.  Geometry  Assumed  in  High  Altitude 
Irirarcd  Radiance  Model 


that  radiates  as  a  blackbody,  and  3)  other  molecules  between  the  lower  and  upi'er 
boundaries.  Radiative  transfer  functions  are  needed  to  compute  the  effects  of 
absorption  of  radiation  by  other  molecules  between  each  of  these  sources  and  the 
molecules  at  point  P.  In  the  absence  of  collisional  excitation  and  deexcitation, 
the  simplest  integral  equation  that  can  be  written  to  describe  radiative  transport 
in  a  finite  plane  parallel  atmosphere  is 


T 

max 

J(T)  =  ^  I  K  j  (|T-t  j )  J(t)  dt  +  S(t)  . 

I' 


(1) 


This  equation  assumes  that  there  is  only  one  radiative  transition  involved  and  the 
radiation  emitted  by  a  molecule  after  it  has  absorbed  radiation  is  independent  of 
the  direction  of  incidence  of  the  absorbed  radiation.  Here  t  is  a  measure  of  the 
optical  path  length  in  the  vertical  direction.  J(t)  is  a  dimensionless  measure  of 
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the  population  of  the  upper  (emitting)  molecular  state.  Kj(t)  is  the  radiative 
transfer  function  that  accounts  for  absorption  of  radiation  emitted  by  molecules 
at  a  distance  t  from  the  absorbing  molecule.  S(t)  is  an  auxiliary  function  that 
includes  the  effects  of  radiation  originating  outside  the  atmosphere.  If  there  is  a 
blackbody  source  located  at  the  lower  boundary  from  which  the  optical  path  length 
is  measured  (r  =  il  at  the  lowei-  boundary),  then 

S(r)  =  1  (t)  .  (2) 

Here  N  is  the  radiance  of  the  lower  boundary  in  arbitrary  units,  such  as  cm 

- 1  ^  j 

ster  ,  or  the  spectral  radiance  in  units  such  as  W  cm  “  ster  (wavelength  inter¬ 
val)  .  K.j(t)  is  the  radiative  transfer  function  that  accounts  for  absorption  be¬ 
tween  the  lower  boundary  and  a  molecule  at  optical  depth  t  above  that  boundary. 

When  a  collimated  beam  of  light  is  incident  on  the  top  of  the  atmosphere  another 
term  must  be  added  to  the  futu  tion  S(t)  in  Kqs.  (1)  and  (2).  If  the  flux  is  incident 
at  a  zenith  angle  whose  cosine  is  u  and  has  an  intensity  in  arbitrary  units,  such  as 
F  W  cm  “  normal  to  its  direction,  the  term  added  to  blr)  has  the  form 
F  M,(t  .  -  T)ftu.  Here  M,(t)  is  the  transmission  along  an  optical  path  with 

i  ITl  (HX  I 

depth  7 . 

The  starting  point  of  the  treatment  of  radiative  transfer  for  the  model  is  the 
treatment  of  a  single  spectral  line.  To  extend  this  treatment  to  a  vibration- 
rotation  band  the  following  assumptions  are  made  (in  addition  to  assuming  that 
each  line  has  the  same  Doppler  profile  shape): 

1.  Hotational  levels  within  a  vibration  level  remain  in  equilibrium  with  the 
translational  temperature.  This  is  valid  up  to  an  altitude  of  at  least  ll(t  km 
according  to  Ooody  (10fi4). 

2.  There  is  complete  redistribution  of  absorbed  energy  among  all  rotational 
levels. 

3.  The  spectral  distribution  of  radiation  from  a  line  is  independent  of  the 
energy  of  the  radiation  absorbed  by  the  level. 

4.  Line  strengths  are  used  as  if  each  line  were  at  the  band  center,  neglecting 
wavelengths  variation  across  a  band. 

.').  There  is  no  overlap  of  adjac  ent  lines. 

(I.  Temperature  variations  at  different  altitudes  do  not  affect  the  radiativi 
transfer. 

These  assumptions  are  introduced  to  make  the  computer  calculation  possible 
in  a  reasonable  length  of  time.  They  do  cause-  tlie  results  to  be  uncertain  to  factors 

fi.  Goody.  R.  M.  (I.'>fi4)  Atmospheric  Radiation.  I,  Theoretical  Basis.  Oxford 
University. 


of  two  and,  thus,  cai’c  must  be  used  in  interpreting  the  results.  Kor  example,  the 
sixth  assumption  results  in  vibrational  populations  of  that  are  too  low  above 
the  mesopause  (80-90  km)  and,  consequently,  computed  radiances  in  the  1.')  pm 
region  may  be  low  by  as  much  as  50  percent. 

In  order  to  simplify  the  computer  programming,  a  single  set  of  radiative  trans¬ 
fer  functions  is  used  for  all  parallel  bands  of  linear  molecules  and  another  set  is 
used  foi'  perpendicular  bands,  rhese  functions  vary  with  the  rotational  temperature. 
In  order  to  account  partially  for  temperature  variation,  the  program  does  not  use 
the  parameter  t,  but  rather  a  derived  parameter  -r^  that  is  the  optical  thickness  at 
the  center  of  the  strongest  line  in  the  band.  This  is  eithei-  an  H  branch  line  for  a 
parallel  band  or  a  Q  branch  line  for  a  perpendiculai-  band. 

In  practice,  rather  than  formally  integrating  a  modified  Kq.  (l),  including  a  t 
dependent  function  to  account  for  coUisional  cleexcitation  and  an  additional  source 
term  to  account  for  coUisional  excitation,  tlu  lOtTiputer  program  sets  an  initial 
vibrational  population  at  each  altitude  and  computt  s  (at  eaih  altitude)  the  total 
radiative  excitation  that  results  fi-om  the  sum  of  radiation  from  thi  earth  and  lower 
atmosphere,  plus  tlu  sun  and  upper  atinosplure.  Vibrational  pop'ulations  are  then 
determined  by  balancing  coUisional  and  radiativ*  xcitation  and  deext  itation  rates. 
Total  band  intensities  are  computed  at  an  oosi  rvation  point  inside  or  outside  the 
atmosphere  by  determining  the  optical  path  from  each  altitude  to  the  observation 
point  and  integrating  the  produc  t  of  loc  al  emission  rale  and  the  band  radiative 
transfer  function  Mj(t). 

The  spectral  distribution  of  radiation  is  determined  by  finding  a  mean  effective 
temperature  and  line-of-sight  column  dinsity  for  the  radiating  molecules,  comput¬ 
ing  intensity  of  each  line  in  the  band  witJi  tlie  liiu  function  I.(t),  summing  the  line 
emission  rates,  and  scaling  each  line  emission  rate  so  that  the  total  band  intensity 
is  obtained. 

Input  data  to  the  computer  program  include  altitude  profiles  of  the  temperature 
and  the  N2,  O.^,  O,  H.^O,  NO,  and  concenti-ations  at  altitudes  between 

ground  level  and  700  km.  N., ,  ©2,  and  t'02  concentrations  below  90  km  are  com¬ 
puted  directly  from  the  temperature  profile.  The  O  concentration  is  based  on  the 
U.  S.  Standard  Atmosphere  1978  (USSA  7G),  but  any  other  source  may  be  used. 

Depending  on  time  and  location  O  concentrations  in  the  critical  region  90  to  12  0  km 

7 

may  vary  by  as  much  as  an  order  of  magnitude  (Offerman  et  al,  1981  );  however, 
this  variation  is  presently  not  included  in  the  program.  Concentrations  of  H2O,  NO. 
and  Og  may  be  uncertain  by  factors  of  two  to  five  or  more  due  to  the  large  expected 
variability  of  these  species  and  the  inaccuracy  of  currently  available  measurements 
(see  Section  3.  5). 

7.  Offerman,  D.,  Friedrich,  V.,  Ross,  P.  and  vonZahn,  V.  (1981)  Neutral  gas  com¬ 
position  measurements  between  80  -  120  km.  Planet  Space  Sci.  (in  publication). 
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There  have  been  relatively  few  measurements  of  the  altitude  variation  of  the 
concentrations  of  infrared  emitting  species  at  altitudes  above  40  km,  the  limit  of 
direct  sampling  by  balloons.  While  ozone  measurements  by  different  scientists 
generally  agree  with  each  other  and  with  theoretical  computations  for  altitudes 

Q 

up  to  100  km  (Weeks  et  al,  1078  )  and  NO  profiles  appear  reasonably  consistent 
^  9 

between  about  80  and  150  km  at  middle  latitudes  (Cravens  and  Stewart,  1978  ), 

corresponding  data  are  less  abundant  for  other  important  species.  There  are  a 
few  cases  where  indirect  data  are  available,  which,  when  combined  with  compu¬ 
tational  models,  can  yield  limits  to  concentrations.  Thus,  Anderson  (197  la)^^ 
has  been  able  to  put  an  upper  limit  on  the  total  OH  abundance  above  60  km  from 
measurement  of  its  resonance  fluorescence.  This  result  allows  a  computation  of 
possible  H2O  profiles  between  50  and  100  km  using  a  photo -chemical  model  that 
includes  vertical  transport  by  eddy  diffusion.  In  a  later  paper  Anderson  (197 Ib)^^ 
gives  a  measured  altitude  profile  for  OH  with  error  limits  of  a  factor  of  two. 

At  present,  estimates  of  chemical  abundances  of  minor  species  important  in 
predicting  infrared  radiation  require  the  use  of  photochemical  computations.  De¬ 
spite  the  inclusion  of  the  best  estimates  of  the  effects  of  transport  processes  such 
as  eddy  mixing,  these  computations  must  be  used  with  some  caution.  The  current 
infrared  model  relies  on  published  measurements  and  recent  calculations  to  esti¬ 
mate  the  abundance  of  infrared  emitting  species  in  the  atmosphere. 

Water  vapor  abundances  are  not  well  known  at  high  altitudes,  even  in  the 
stratosphere.  Direct  measurements  of  water  vapor  abundances  place  lower  and 
upper  limits  of  the  volume  mixing  ratio  at  two  and  ten  ppm  in  the  upper  strato¬ 
sphere.  In  a  recent  optical  measurements  from  a  balloon  Goldman  et  al  (1973) 
observed  the  thermal  emission  from  the  rotational  spectrum  near  25  pm  and 
determined  that  the  volume  mixing  ratio  near  30  km  was  nearly  three  ppm. 


8.  Weeks,  L.  H.,  Good,  R.E.,  Randhawa,  J.S.  and  Trinks,  H.  (1978)  Ozone 

measurements  in  the  stratosphere,  mesosphere,  and  lower  thermosphere 
during  Alladdin  74,  J.  Geophys.  Res.  83.-978-982. 

9.  Cravens,  T,  E.  and  Stewart,  A.  I.  (1978)  Global  morphology  of  nitric  oxide 

in  the  lower  E  region,  J.  Geophys.  Res.  83:2446-2452. 

10.  Anderson,  J.G.  (197  la)  Rocket-borne  ultraviolet  spectrometer  measurement 

of  OH  resonance  fluorescence  with  a  diffusive  transport  model  for  meso¬ 
spheric  photochemistry,  J.  Geophys.  Res.  76:4634-4652. 

11.  Anderson,  J.G.  (1971b)  Rocket  measurement  of  OH  in  the  mesosphere, 

J.  Geophys.  Res.  7j6:7820-7824. 

12.  Goldman,  A,,  Murray,  D.  G. ,  Murcray,  F,  H.,  Williams,  W.J.  and 

Brooks,  J.N.  (1973)  Distribution  of  water  vapor  in  the  stratosphere  as 
determined  from  balloon  measurements  of  atmospheric  emission  spectra 
in  the  24-29  pm  region,  Appl.  Opt.  12:1045-1053. 
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Radford  et  al  (1977)  used  ground  based  microwave  soundings  to  determine  meso¬ 
spheric  water  vapor  concentrations.  They  obtained  10  ppmV  at  50  km  v.  ith  a  max¬ 
imum  of  15  ppmV  at  60  km.  Rogers  et  al  (1977a)^^  have  analyzed  radiance  meas¬ 
urements  obtained  from  the  AFGL/DNA  ICECAP  program  to  infer  an  H2O  mixing 
ratio  of  3.  5  ppmV  between  50  and  70  km  during  late  winter  at  Fairbanks,  Alaska. 

At  present,  5  ppmV  is  assumed  as  a  reasonable  mixing  ratio  for  use  in  the  radi¬ 
ance  model.  Above  90  km  this  value  is  too  high  because  of  photodissociation.  Dur¬ 
ing  summer  at  high  latitudes,  a  value  one  tenth  of  this  may  be  appropriate  near  the 
cold  mesopause.  Figure  2  compares  the  model  calculations  of  the  H2O  radiance  at 

6.  3  pm  with  the  AFGL.  SPIRE  measurements  from  a  rocket  payload  launched  at 

5 

Poker  Flat  on  28  September  1977  (Stair  et  al,  1981  ).  There  is  good  agreement 
between  60  and  70  km  where  the  model  and  data  overlap.  However,  the  AFGL  ELE 
experiment,  launched  from  White  Sands  on  6  August  1977  requires  that  the  .^2^ 
concentrations  be  reduced  by  a  factor  of  five  at  90  km  in  order  to  achieve  agreement. 
These  differences  may  simply  represent  daily,  seasonal,  or  latitudinal  variations 
in  the  atmospheric  water  vapor  content. 

A  large  number  of  measurements  of  NO  concentrations  are  available.  Repre¬ 
sentative  midlatitude  concentrations  are  shown  in  Figure  3  that  are  based  on  the 

15 

twilight  measurements  of  Baker  et  al  (1977)  and  the  calculations  of  Oran  et  al 
1 6 

(1975)  ,  NO  concentrations  are  greater  at  higher  latitudes.  The  nighttime 

SPIRE  data  at  5.4  pm  require  the  use  of  the  higher  daytime  concentrations  of 
Figure  3,  indicating  that  the  production  of  NO  is  enhanced  at  higher  latitudes 
(see  Section  3.  5). 

Calculations  of  Og  abundances  must  take  into  account  lateral  transport  bv 
winds,  as  well  as  vertical  transport  by  eddy  mixing  and  molecular  diffusion,  in 
addition  to  photochemistry.  Current  photochemical  models  agree  with  observa¬ 
tions  obtained  by  measuring  absorption  of  solar  and  stellar  ultraviolet  radiation 
and  those  deduced  from  measurements  of  the  molecular  oxygen  1.  37  uin  (O., 
infrared  atmospheric  band.  The  exact  density  profile  depends  strongly  on  assumed 
eddy-diffusion  coefficients  because  of  the  chemistry  involving  atomic  H  and  O 

13.  Radford,  H.E.,  Litvak,  M.M.,  Gottlieb,  C.A.  Rosenthal,  S.  K.  and  Lilley, 

A.E.  (1977)  Mesospheric  water  vapor  measurements,  J.  Geophys.  Res. 
^:472-478. 

14.  Rogers,  J.W.,  Stair,  A.T. ,  Degges,  T.  C. ,  Wyatt,  C.  L.  and  Baker,  U.  J. 

(1977)  Rocketborne  measurement  of  mesospheric  H^O  in  the  auroral  zone, 
Geophys.  Res.  Lett.  4:366-368. 

13.  Baker,  K.  D. ,  Nagy,  A.  F. ,  Olsen,  R.O. ,  Oran,  E.S.,  Randhawa,  J., 

Strobel,  D.  F.  and  Tohmatsu,  T.  (1977)  Measurement  of  the  nitric  oxide 
altitude  distribution  in  the  mid-latitude  mesosphere,  J.  Geophys.  Res. 
^:3281-3286. 

16.  Oran,  J.S. ,  Julieene,  P.S.  and  Strobel.  D.  F.  (1975)  Odd  nitrogen  in  the 
thermosphere,  J.  Geophys.  Res.  80:3068-3076. 


Figure  2.  Comparison  of  SPIRE  Water  Vapor  Measurement  at 
6.  5  nm  With  Model  Calculations 


concentrations.  Og  abundances  used  in  the  current  infrared  radiance  model  are 

shown  in  Figure  4.  Below  GO  km  the  ozone  model  in  the  USSA  76  is  used.  These 

17 

profiles  rely  heavily  on  measurements  by  Llewellyn  and  Wilt  (1977)  and  Hays 
and  Roble  (1073).  The  calculations  of  Thomas  and  Bowman  (1972)^^  have  been 
used  as  a  guide  to  predict  diurnal  variations  in  the  O.^  abundances.  Diffusive 
equilibrium  is  assumed  above  100  km. 

Figure  5  compares  SPIRE  measurement  data  for  9.  6  pm  Og  limb  viewing 
radiances  and  those  computed  with  the  Og  profiles  of  Figure  4.  Such  good  agree¬ 
ment  is  to  some  extent  fortuitous.  In  matching  Og  radiance  profiles  taken  from 
zenith  viewing  rocketborne  instruments  with  model  predictions,  features  of  the 
nighttime  Og  concentration  profile  must  be  shifted  at  times  by  as  much  as  5  km  in 
altitude.  Concentrations  may  vary  by  as  much  as  a  factor  of  three. 

17.  Llewellyn,  E,  J,  and  Witt,  G.  (1977)  The  measurement  of  ozone  concentrations 

at  high  latitude  during  the  twilight.  Planet.  Space  Sci.  25: 165-172. 

18.  Hays,  P.B.  and  Roble,  R.G.  (1973)  Observations  of  mesospheric  ozone  at 

low  latitudes.  Planet  Space  Sci.  21:273-279. 

19.  Thomas,  L.  and  Bowman,  M.R.  (1972)  The  diurnal  variations  of  hydrogen 

and  oxygen  constituents  in  the  mesosphere  and  lower  thermosphere, 

J.  Atmos.  Terr.  Phys.  34: 1843-1858. 
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ALTITUDE,  km 


Figure  4,  O3  Concentration  Profiles  Used  in  the  High  Al 
titude  Infrared  Radiance  Model 


Figure  5.  Comparison  of  SPIRE  9.6  pm  Data  With  Model  Calculations 

COg  is  dissociated  by  solar  radiation  and  lost  by  reaction  with  in  the  upper 
atmosphere.  Below  about  110  km  it  is  well  shielded  from  dissociation  by  molecu¬ 
lar  oxygen  and,  thus,  should  remain  well  mixed  below  100  km.  COg  concentrations 
have  been  reported  for  two  sets  of  measurements  using  mass  spectrometers. 
Philbrick  et  al  (1973)  presented  measurements  that  showed  that  the  CO^  was  in 
mixing  equilibrium  between  85  and  100  km.  The  measurements  by  Trinks  and 
Fricke  (1978)  indicated  that,  if  the  nighttime  profile  is  assumed  to  be  in  dif¬ 
fusive  equilibrium  above  100  km,  a  calculation  that  accounts  for  dissociation  and 
loss  of  CO2  by  photodissociation  and  reaction  with  O  gives  reasonable  agreement 
for  the  daytime  lower  stratosphere.  The  current  model  does  not  include  the  day- 
time  loss  processes  implied  by  the  results  of  Trinks  and  Fricke  (1978).  Com¬ 
plete  mixing  is  assumed  up  to  a  level  between  90  and  120  km  and  diffusive  equilib¬ 
rium  is  assumed  atwve  that  variable  level. 


20.  Philbrick,  C.R.,  Faucher,  G.A.  and  Trzcinski,  E.  (1973)  Rocket  measure¬ 

ments  of  mesospheric  and  lower  thermospheric  composition.  Space  Research 
Xm;255-260.  - 

21.  Trinks,  H.  and  Fricke,  K.  H,  (1978)  Carbon  dioxide  concentrations  in  the 

lower  thermosphere,  J.  Geophys.  Res.  83:3883-3886. 
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Figure  6  compares  the  measured  SPIRE  15  pm  CO2  radiance  and  the  model 
calculations.  Agreement  is  good  below  95  km,  but  above  that  altitude  the  meas¬ 
ured  radiance  levels  are  about  twice  the  model  predictions.  Part  of  the  discrep- 
ancy  is  due  to  use  of  assumption  6  above.  Kumer  (1979)  has  obtained  better 
agreement  by  treating  the  temperature  variation  near  the  mesopause,  but  that 
analysis  still  finds  the  need  for  another  excitation  mechanism.  Sharma  and  Nadile 
(1980)  have  also  analyzed  this  data.  They  obtai-.ed  agreement  between  the  data 

and  predictions  by  assuming  an  oxygen  concentration  with  a  peak  value  of 
11  -3 

1.7X10  cm  at  105  km  and  a  room  temperature  deactivation  rate  of  CO„  by 

•  13  3  ^ 

oxygen  atoms  of  5  X  10  cm  /sec. 

Figure  7  shows  a  comparison  of  :ho  predictions  of  the  High  Altitude  Radiance 
Model  and  the  SPIRE  data  over  tht  total  spectral  range  measured.  In  the  spectral 
region  from  9  to  17  pm  the  predi  '.ums  ap.,ear  to  be  too  iow. 


Figure  6.  Comparisonof  SPIRE  15pm  CO,  Radiances  With  Model 
Calculations 


22.  Kumer,  J.  B.  (1979)  Private  communication. 

23.  Sharma,  R.D.  and  Nadile  R.M.  (1981)  Carbon  dioxide  (1^2)  radiance  results 

using  a  new  non-equilibrium  model.  Paper  presented  at  AIAA  Aerospace 
Sciences  Meeting.  Jan.  1981,  St.  Louis.  Missouri,  Paper  no. 
AIAA-81-0421. 


WAVELENGTH  (um) 


Figure  7.  Comparison  of  SPIRE  Measurements  at  90  km  of  the  Sunlit  Exoatmospher- 
ic  LimbViewing  RadianceWith  Model  Calculations  Under  Sunlit  Conditions  at  100  km 

The  largest  uncertainties  in  applying  the  infrared  radiance  model  result  from 
lack  of  knowledgeof  the  atmospheric  temperature  and  molecular  species  concentra¬ 
tions.  Further  uncertainties  arise  from  the  fact  that  not  all  radiating  molecules 
or  possible  excitation  mechanisms  are  included  in  the  program.  Also,  other  bands 
of  those  molecules  that  are  included  may  have  significant  effects.  The  modeling 
effort  relies  on  experiments  to  identify  the  most  appropriate  additions  for  future 
work. 
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I.  IMTROUliC.TION 

As  discussed  In  Chapter  2,  predictions  of  the  infrared  radiance  are  made 
principally  on  the  basis  of  theoretical  models  that  utilize  observations  of  the 
atmospheric  temperatures  and  species  concentrations.  Extensive  measurements 
of  these  properties  of  atmosphere  have  been  carried  out  over  the  past  ten  years 
using  ground-based,  rocket  and  satellite  borne  instruments.  F'rom  these  measure¬ 
ments  a  number  of  models  have  been  developed  to  describe  the  behavior  of  the 
smooth  atmospheric  profiles  that  result  from  variations  in  the  solar  heating,  mag¬ 
netic  activity,  and  motion  of  the  earth. 

In  this  chapter  these  models  and  measurements  that  support  these  models  are 
reviewed  briefly.  Special  emphasis  is  given  to  the  Jacchia  77  (J77)  model  as  this 
model  has  been  selected  as  the  basis  for  updating  the  high  altitude  infrared  model 
to  include  the  effects  of  atmospheric  variability.  This  review  will  enable  us  to 
draw  some  inferences  concerning  the  variability  of  the  atmospheric  infrared 
radiance  that  can  result  from  the  variations  in  species  concentrations  and 


temperature.  In  Chapter  4  initial  efforts  to  include  some  of  the  variations  repre¬ 
sented  by  these  models  in  the  calculation  of  the  infrared  radiance  arc  presented. 


2.  ATMOSPHERIC  MODELS 

2.1  C1RA72 

In  CIRA  72  (COSPAR,  1972  ^  the  latest  internationally  adopted  model,  a  mean 
COSPAR  International  Reference  Atmosphere  compiled  by  K.  S,  Champion  and 
R..A.  Schweinfurth  for  the  altitude  range  25  to  .500  km  is  presented.  This  atmos¬ 
phere  represents  average  diurnal,  seasonal,  semiannual,  and  geomagnetic  con- 

-  .  9  -  o 

ditions  appropriate  to  30°N  latitude  and  a  10.7  cm  solar  flux  of  145  ^  lO  ““W  m 
Hz  The  resultant  exospheric  temperature  is  1000  K.  In  addition,  CIRA  72  pre¬ 
sents  tables  of  the  atmospheric  properties  as  functions  of  latitude,  time  of  year 
for  altitudes  from  25  to  110  km  due  to  G.V.  Groves  and  tables  of  atmospheric 
profiles  for  altitudes  from  110  to  1000  km  with  exospheric  temperatures  ranging 
from  500  to  2200  K,  derived  from  the  work  of  L,  Jacchia.  Experimental  data 
that  support  these  models  are  also  discussed,  as  well  as  the  variations  to  be  ex¬ 
pected  in  these  models. 

2.2  IJSSA76 

The  U.  S.  Standard  Atmosphere  1976,  (USSA  76)  (COESA,  1976^)  presents  a 
single  consistent  model  of  globally  averaged  values  of  the  atmospheric  parameters 
for  altitudes  up  to  1000  km.  This  model  is  assumed  to  be  representative  of  moder¬ 
ate  solar  activity.  The  presentation  of  this  model  included  a  review  of  some  fea¬ 
tures  of  the  variation  of  atmospheric  properties.  Seasonal  and  latitudinal  vari¬ 
ations  in  the  mean  monthly  temperature  and  density  below  90  km  are  compared 
with  average  values  of  USSA  76  in  Figures  1  and  2  and  are  in  general  agreement 
with  similar  results  presented  in  CIRA  72.  Figure  2  shows  that  variations  of  as 
much  as  50  percent  in  density  can  occur  at  altitudes  at  which  significant  infrared 
radiation  is  produced.  Such  variations  in  temperature  and  density  can  lead  to 
appreciable  variations  in  the  infrared  radiance.  The  variation  in  species  con¬ 
centrations  which  result  from  variations  in  the  temperature-height  profiles  at 


1.  COSPAR  (1972)  Committee  on  Space  Research,  CIRA  72,  Cospar  International 

Reference  Atmosphere.  Akademle-Verlag,  Berlin. 

2.  Jacchia,  L.G.  (1977)  Thermospheric  Temperature,  Density  and  Composition: 

New  Models,  SAO  Special  Report  No.  '675, 
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altitudes  above  100  kin  for  various  degrees  of  solar  and  geomagnetic  activity  are 
also  discussed  in  USSA  76. 


2.3  Air  Kori-e  Reference  Alinosplieres 

Recently  Cole  and  Kantor  (1978)  published  reference  atmospheres  that  pro¬ 
vide  seasonal,  latitudinal,  and  day-to-day  variability  of  the  thermodynamic  pro¬ 
perties  of  the  atmosphere  up  to  an  altitude  of  90  km.  These  models  expand  upon 
and  update  information  contained  in  the  USSA  76  (COESA,  1976^).  Cole  and 

4 

Kantor  present  sets  of  mean  monthly  reference  atmospheres  describing  seasonal 
changes  in  the  vertical  distributions  of  temperature,  density,  and  pressure  for 
15°  intervals  of  latitude  between  the  equator  and  pole.  In  addition,  they  present 
estimates  of  the  day-to-day  variations  around  the  monthly  median  values  of  tem¬ 
perature  and  density.  Figure  3  displays  the  latitudinal  variations  in  average 
monthly  temperature  for  the  months  of  January  and  July.  From  these  temperature 
profiles  density  profiles  can  be  constructed  using  the  barometric  equation.  The 
results  of  this  calculation  show  density  variations  similar  to  those  shown  in 
F igure  2 . 

In  arctic  and  subarctic  regions,  sudden  warmings  and  coolings  of  the  winter 
stratosphere  and  mesosphere  produce  large  changes  in  the  vertical  and  horizontal 
structure  of  the  atmosphere,  both  in  the  magnitude  and  in  the  altitude  region  of 
the  temperature  and  density  fluctuations.  The  di._  -to-day  variations  in  temper¬ 
ature  are  in  some  cases  as  great  or  greater  than  the  seasonal  or  latitudinal 
changes.  Although  warmings  and  coolings  occur  throughout  the  artic  and  sub- 
artic  regions,  the  largest  changes  generally  occur  between  latitudes  60°  and  7  5°. 

2.4  MSIS 

The  MSIS  model  (Hedin  et  al,  1977a,  ^  1977b*’)  is  a  global  thermospheric 
model  that  has  been  constructed  from  measurements  of  the  density  of  N2  with 
mass  spectrometers  on  five  satellites  (AE-B,  OGO  6,  San  Marco  3,  Aeros  A,  and 
AE-C)  and  from  neutral  temperatures  inferred  from  incoherent  radar  backscatter 
measurements  at  four  ground  stations  (Arecibo,  Jicamarca,  Millstone  Hill,  and 

4.  Cole,  A.E.  and  Kantor,  A.  J.  (1978)  Air  Force  Reference  Atmospheres, 

AFGL-TR-78-0051,  AD  058  505. 

5.  Hedin,  A.E.,  Salah,  J.E.,  Evans,  J.V.,  Reber,  C.A.  ,  Newton,  G.  P. , 

Spencer,  N.W.,  Kayser,  D.C.,  Alcayde,  D. ,  Bauer,  P.  ,  Cogger,  L. ,  and 
McClure,  J.P.  (1977a)  A  global  thermospheric  model  based  on  mass  spec¬ 
trometer  and  incoherent  scatter  data,  MSIS  1,  N2  density  and  temperature, 

J.  Geophys.  Res.  82:2139-2147. 

6.  Hedin,  A.E.,  Reber,  C.A.,  Newton,  C.  P. ,  Spencer,  N.W. ,  Brinton,  H.  C. , 

and  Mayr,  H.G.  (1977b)  A  global  toermospheric  model  based  on  mass  spec¬ 
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82:2148-2156. 


Figure  3.  Latitudinal  Temperature-height  L'ross  Sections  of  Monthly  Temperatures  for  January  and  July  <C'ole  and 
Kantor,  1978'*) 


St.  Santin).  The  satellite  and  ground-based  incoherent  scatter  measurements 
provide  unique  and  complementary  information  on  the  structure  of  the  thermos¬ 
phere.  In  addition,  the  MSIS  model  combines  measurements  of  the  densities  of 
O,  He,  and  Ai  from  neutral-gas  mass  spectrometers  on  four  satellites  and  in¬ 
ferred  ©2  and  H  densities  from  an  ion  mass  spectrometer  on  the  AE-C  satellite 
with  the  neutral  temperature  and  N2  density  models  to  produce  a  global  model  of 
thermospheric  composition  in  terms  of  inferred  variations  at  120  km.  The  over¬ 
all  data  set  of  these  measurements  covers  the  time  period  from  the  end  of  1005  to 
mid- 197 5.  In  MSIS  the  exospheric  temperature  and  other  quantities  are  repre¬ 
sented  by  least  squares  fits  to  the  data  in  spherical  harmonic  expansions.  In  197  9 
the  MSIS  model  was  extended  to  include  longitudinal  and  universal  time  (UT)  vari- 
ations  in  thermospheric  neutral  temperature  and  composition  (Hedin  et  al,  1979  ). 
Terms  that  depend  only  on  longitude  indicate  a  temperature  enhancement  at  12  0  km 
of  about  30°K  at  a  latitude  of  approximately  70°  along  the  longitude  of  the  magnetic 
poles.  As  with  the  variations  with  magnetic  activity,  the  extrapolated  120  km  vari¬ 
ations  are  generally  in  phase  with  temperature  variations  for  N2,  Ar,  and  and 
out-of-phase  for  He,  O,  and  H.  The  co.mbined  longitude  and  UT  variations  reflect 
the  influence  of  the  earth's  magnetic  field  and  suggest  that  the  primary  variations 
may  be  simply  represented  in  magnetic  coordinates. 

The  MSIS  model  is  most  representative  of  magnetically  qniet  conditions  with 
the  best  data  coverage  between  altitudes  200  and  600  km  and  with  solar  ^ 
indices  between  75  and  180.  The  species  concentrations  should  be  reasonably 
valid  above  150  km,  but  below  150  km  they  are  essentially  extrapolations  from 
higher  altitudes.  Temperatures  below  200  km  may  not  adequately  represent  di¬ 
urnal  and  semidiurnal  variations.  The  limitations  on  altitudes  make  the  applica¬ 
tion  of  the  MSIS  model  to  the  calculation  of  limb  viewing  infrared  calculations  very 
difficult.  In  addition,  the  use  of  spherical  harmonic  construction  of  this  model 
make  it  difficult  to  join  to  the  lower  altitude  atmospheric  models. 

2.5  Jarchia  77  (J77) 

2 

The  Smithsonian  Astronomical  Observatory  models  (Jacchia,  1977  )  that 
describe  the  thermospheric  temperature,  density,  and  composition  consist  of 
essentially  two  parts;  1)  the  basic  static  models,  that  give  temperature  and  den¬ 
sity  profiles  of  the  atmospheric  constituents  for  any  specified  exospheric  temper¬ 
ature  and  2)  a  set  of  formulae  for  correcting  the  exospheric  temperature  based  on 
various  types  of  thermospheric  variation. 


7.  Hedin,  A.E.,  Reber,  C.  A.  ,  Spencer,  N.W.,  Brinton,  H.C'.  and  Kayser,  D.  C. 
(1979)  Global  model  of  longitude/UT  variations  in  thermospheric  composition 
and  temperature  based  on  mass  spectrometer  data,  J.  Geophys.  Res. 
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The  Jacchia  models  are  based  primarily  on  a  very  large  number  of  satellite 
drag  measurements  and  on  composition  measurements  from  several  satellites, 
OGO-6,  ESRO  4,  and  S3-1.  For  the  basic  static  models,  tables  of  temperature, 
species  concentrations,  pressure  and  mass  density  are  given  for  altitudes  from 
90  to  2500  km  and  for  exospheric  temperatures  from  500  to  2600  K.  The  tem¬ 
perature  profiles  start  from  a  constant  value  of  188  K  at  an  altitude  of  90  km, 
rise  to  an  inflection  point  of  125  km,  and  become  asymptotic  to  an  "exospheric" 
temperature  at  approximately  500  km.  Examples  of  the  resultant  temperature 
profiles  are  shown  in  Figure  4.  The  mass  density  is  assumed  to  have  a  constant 
value  of  3.43  X  10  kg/m  at  the  altitude  of  90  km.  The  atmosphere  is  assumed 
to  be  composed  of  Ng,  O2,  O,  Ar,  He,  and  H  with  constant  mixing  ratios  up  to 
100  km  and  with  abundances  determined  by  diffusive  processes  above  this  height. 
Number  density  profiles  of  the  major  atmospheric  species  for  two  exosphere 
temperatures  are  illustrated  in  Figure  5.  The  assumption  of  fixed  boundary  con¬ 
ditions  at  90  km  can  lead  to  discontinuities  in  the  density  and  temperature  profiles 
if  the  Jacchia  model  is  joined  to  atmospheric  models  that  describe  the  mesosphere. 
Indications  of  the  possible  density  variability  and  temperature  variability  at  90  km 
"■re  shown  in  Figures  1,  2,  and  3. 


Figure  4.  Foiu"  Temperature  Profiles  from  the  J77  Model 
(Jacchia,  19772) 
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Figure  5.  Variations  of  Atmospheric  Species  Concentrations  at  78°  Latitude 
as  Predicted  by  J77  Model  Under  Two  Conditions  of  Geomagnetic  Activity 
and  Exospheric  Temperature 


The  types  of  variation  that  the  J77  models  attempt  to  represent  can  be  classi¬ 
fied  as  follows: 

1.  Variation  with  solar  cycle 

2.  Variation  with  changes  in  activity  on  the  disk  of  the  sun 

3.  The  daily  or  diurnal  variation 

4.  Variations  with  geomagnetic  activity 

5.  Seasonal-latitudinal  variations 

6.  The  semiannual  variation 

All  these  variations  are  based  on  parameters  that  can  be  predicted  with  vary¬ 
ing  degrees  of  accuracy  or  can  be  obtained  from  observational  data.  Such  static 
models  should  be  quite  adequate  to  describe  the  mean  conditions  when  the  charac¬ 
teristic  time  of  the  variation  is  much  longer  than  the  time  scale  of  conduction, 
convection,  and  diffusion  processes.  The  following  paragraphs  describe  how  the 
above  variations  are  incorporated  by  Jacchia. 

The  magnitude  of  the  variation  of  the  exospheric  temperature  is  represented 
as  a  function  of  both  the  instantaneous  value  of  the  10.  7 -cm  solar  flux  F  and  the 
flux  averaged  over  a  few  solar  rotations  F.  The  variation  of  exospheric 
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tfnifK  i'aturc  is  applied  with  a  time  lag  that  varies  from  0.  f>  days  at  noon  local 
solar  time  (LST)  to  l.(i  days  at  midnight  l^ST. 

The  diurnal  variation  of  tl\e  exosplierie  temperature  is  represented  as  the 
sum  of  two  terms,  one  repiesenting  a  seasonal-latitudinal  variation  that  is  in  - 
pendent  of  local  time  and  another  that  is  a  function  of  local  time.  The  seasonal- 
latitudinal  (or  annual)  variation  has  a  maximum  value  of  1')  percent  at  tlie  poles 
and  is  proportional  to  tlic  declination  of  the  sun  and  to  the  sine  of  the  latitude.  In 
Tigure  r.  the  resulting  variations  of  exospheric  lemperaUire  at  the  equinoxes  and 
the  June  solstice  an  shown  for  quiet  magnetic  conditions  =  0).  This  global 
variation  of  exospheric  temperature,  in  conjunction  with  the  resulting  variatv':. 
in  species  concentrations,  can  produce  a  variation  in  the  global  infrared  emission. 
Examples  of  such  variations  are  shown  i.i  Chapter  -1. 

The  variation  with  geomagnetic  activity  is  repi'oscntcd  as  a  temperatuie 
increase  that  is  a  function  of  the  geomagnetic  latitude  and  the  index.  Hie 
accompanying  variation  in  species  concentrations  is  comprised  of  three  separate 
components:  n  thermal  component  due  to  the  change  in  scale  height  induced  by 
the  ihange  in  temperature,  a  component  due  to  the  change  in  the  altitude  of  the 
homopmise  (the  interface  between  the  regimes  of  mixing  and  diffusion),  and  a  iom- 
pon>  nt  due  to  the  "equatorial  wave  '  that  results  from  a  density  enhancement  re¬ 
sulting  fi'om  convei'tion  toward  the  equator.  An  example  of  the  resulting  change  in 
exospheric  temperature  is  shown  in  Figure  7,  where  ii  is  seen  that  the  tempei-ature 
at  the  higher  latitudes  can  be  significantly  changed  by  magnetic  activity.  Tlic  ef¬ 
fective  change  in  thi  altitude  of  the  homopause  included  in  the  model  is  shown  in 
Figure  8  as  a  function  of  for  various  geomagnetii’  latitudes.  Hie  density  vari¬ 
ations  of  four  atmospheric  constituents  as  a  function  of  the  invariant  magnetic 
latitude  for  a  =  8  are  shown  in  Figure  fi.  The  effect  of  the  geomagnetic  correc¬ 
tion  to  the  J77  model  is  also  shown  in  tlie  profiles  presented  in  Figure  5. 

A  recent  modification  of  the  J77  model  (Phiibrick  and  Gardner,  1980^^)  pre¬ 
dicts  the  nighttime  altitude  profile  of  C'O,^  as  a  function  of  latitude  and  geomagnetic 
activity.  In  Figure  10  the  variation  in  this  important  infrared  contributor  is  shown 
for  low  and  high  latitudes  and  for  a  range  of  geomagnetic  activity.  These  results 
suggest  that  particularly  significant  variations  in  the  infrared  radiation  from  GO^ 
may  occur  with  the  variation  in  geomagnetic  activity.  Calculations  of  the  radiation 
from  GO.,  presented  in  Chapter  4  demonstrate  this  effect. 


8.  Phiibrick,  G.R.  and  Gardner,  M.E.  (1980)  Private  communication. 


Figure  6.  Global  Distribution  of  the  Exospheric  Temperature  for  Quiet  Geomag¬ 
netic  Conditions  (Kp  =  0),  The  coordinates  are  local  solar  time  and  geographic 
latitude  (Jacchia 
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Figure  7.  Exospheric  Temperature  Profiles  Along  the 
Li  -  rilete  (SfiOO)  Meridional  Circle  Along  Which  tlie  Local 
Solar  Time  is  17"  in  One  Hemisphere  and  5"  in  the  Other, 
for  Various  Levels  of  Ceoma«ietic  Activity.  Even  a  mod¬ 
erate  level  of  activity  (Kp  »  2)  has  the  effect  of  shifting  the 
temperature  maximum  from  the  ^quator  to  the  poles  at  the 
time  of  equinoxes  (Jaeehia,  1!‘77-) 
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3.  ATMOSPHERIC  MEASUREMENTS  COMPARED  WITH  MODELS 

Various  measurements  have  been  compared  to  the  predictions  of  the  above 

models.  These  comparisons  show  that  the  J77  and  MSIS  models  are  in  general 

agreement.  However,  neither  provides  a  completely  satisfactory  representation 

of  the  atmosphere.  For  example,  in  Figure  11  altitude  profiles  of  N„,  Ar,  O,  and 

^  9 

©2  number  densities  measured  during  the  Aladdin  74  program  (Trinks  et  al,  1978' ) 
are  compared  to  the  predictions  of  the  J77  and  MSIS  models.  The  dip  in  the  den¬ 
sities  by  a  factor  of  3  to  5  at  112.  5  km  is  indicative  of  the  variations  in  the  lowest 

p 

part  of  the  thermosphere.  Trinks  et  al  suggested  that  this  variation  could  be  the 
superposition  of  a  semidiurnal  tidal  wave  and  a  gravity  wave.  This  proposition 
was  supported  by  the  temperature -height  profile  obtained  from  incoherent-scatter 
radar  measurements. 

Figures  12  and  13  show  the  comparison  of  model  prediction  and  measurements 
at  higher  altitudes.  Predictions  of  the  N2,  O,  and  Ar  density  models  using  the 
J77  and  MSIS  are  compared  to  measurements  made  with  the  S3-1  satellite 

9.  Trinks,  H.  and  Fricke,  K.H.  (1978)  Carbon  dioxide  concentrations  in  the 
lower  thermosphere,  J.  Geophys.  Res.  83:3883-3886. 
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Figure  11.  Density-height  Profiles  of  N2,  Ar,  O,  and  O2.  The 
dots  and  heavy  solid  curves  represent  the  descent  data;  the  ascent 
data  are  given  by  the  dashed  curves.  The  number  densities  are  in 
units  per  cm^  and  on  a  logarithmic  scale.  Also  shown  are  the 
predictions  of  the  Jacchia  77  model  (J77)  and  the  MSIS  model 
(Hedin,  et  al,  1977a,  ^  1977b®);  the  triangles  denote  rocleet  aver¬ 
ages  (Offerman,  1974I®),  (Trinlis  et  al,  1978®) 


Q 

(Philbrick  and  Gardner  1980  )  for  both  magnetically  quiet  and  disturbed  conditions. 
The  differences  between  the  models  and  the  data  show  that,  while  general  agree¬ 
ment  exists,  the  observed  small-scale  variations  of  the  densities  are  not  modeled. 
It  is  also  evident  from  these  two  typical  cases  that,  during  a  period  of  higher  geo¬ 
magnetic  activity,  such  as  shown  in  Figure  13.  significantly  more  structure  is 
observed  in  the  profiles. 


10.  Offerman,  D.  (1974)  Composition  variations  in  the  lower  thermosphere, 
J.  Geophys.  Res.  79:4281-4293. 
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4.  (;eomagnktic  induced  variations  of  species  concentrations 


The  influence  of  geomagnetic  activity  on  the  upper  atinospherit  species  den¬ 
sities  is  even  stronger  than  it  is  upon  the  total  mass  density.  This  effect  is  due 
to  the  difference  in  diffusion  velocities  of  the  various  species  and  results  in  de¬ 
pletions  of  those  species  with  mass  lighter  than  the  mean  molecuiar  weight  (see, 
for  example,  the  predictions  of  the  Jacchia  77  model  in  Figures  5  and  10).  Fig¬ 
ure  13  shows  that  disturbed  magnetic  conditions  (K^  =  6  )  can  greatly  increase 
N,  density,  Ar  density,  and  the  structure  in  the  density  profiles. 

in  Figure  14  the  correlation  of  species  concentrations  with  geomagnetic  acti¬ 
vity  is  studied  further.  It  is  seen  that  the  Ar/N^  ratio,  as  well  as  the  total  den¬ 
sity,  is  strongly  correlated  with  the  level  of  geomagnetic  activity.  The  Ar  den¬ 
sity  is  highly  correlated  with  the  level  of  geomagnetic  activity  whereas,  the  O 
density  appears  to  have  a  weak  negative  correlation  with  K^,  about  30  percent 
change  in  this  case.  The  large  increase  in  Ar  density  with  Kp  is  shown  further  in 
Figure  15,  where  measurements  made  by  the  S3-1  satellite  are  plotted  versus  geo¬ 
magnetic  latitude  for  two  ranges  of  K^,  0  £  £  2*  and  4  £  K  £  The  depend¬ 

ence  on  geomagnetic  activity  of  CO^  is  expected  to  be  similar  to  that  of  Ar  due  to 
the  nearly  equal  atomic  masses.  Thus,  the  effect  of  increased  geomagnetic  acti¬ 
vity  upon  atmospheric  infrared  radiance  should  be  very  important.  The  influence 
of  photochemical  processes  on  this  connection  is  considered  in  the  next  section. 

In  Figure  16  estimates  of  the  frequency  of  occurence  of  high  magnetic  activity 
and  of  the  duration  of  such  increases  are  shown.  These  estimates  are  from  the 
work  of  Pazick  (1976).  He  presented  tables  of  probabilities  of  geomagnetic 
activity.  In  particular,  these  results  indicate  that  approximately  one-sixth  of  the 
time  Kp  will  exceed  4  and  approximately  one-half  these  disturbances  will  last 
longer  than  one  day. 


I 


5.  VARIATIONS  IN  CONCENTRATIONS  OF  THE  MINOR  ATMOSPHERIC  SPECIES 

I 

Direct  measurements  of  the  minor  atmospheric  species  have  been  carried 
out  from  rocketborne  platforms.  The  results  of  these  measurements  make  it 
possible  to  draw  some  conclusions  concerning  the  overall  variability  of  the  minor 
species,  but  few  conclusions  can  be  drawn  concerning  the  time  dependence  of  these 
variations. 


11.  Pazick,  P.  M.  (1976)  Conditional  Probabilities  of  the  Geomagnetic  Index  A  , 
AFGRL-TR -76-0034,  AD  A023  675.  P 
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Figure  15.  Ar  Density  at  an  Altitude  of  160  km  as  a 
Function  of  Geomagnetic  Latitude  for  Quiet  and  Dis¬ 
turbed  Magnetic  Conditions 


Measurements  of  the  density  of  CO  by  mass  spectrometer  (Philbrick  et  al, 

12  ^ 

1973  )  have  shown  that  COg  is  mixed  in  its  ground  level  ratio  between  85  and 

100  km.  Measurements  of  the  CO„  density  above  the  turbopause  were  performed 

^  13 

at  night  in  1973  by  Offermann  and  Grossman  (1973)  and  during  the  day  in  1974 

12.  Philbrick,  C.R.,  Faucher,  G.  A.  and  Trzcinski,  E.  ( 1973)  Rocket  measure¬ 

ments  of  mesospheric  and  lower  thermospheric  composition.  Space  Research 
Xin:255-260.  - 

13.  Offerman,  D.  and  Grossman,  K.U.  (1973)  Thermospheric  daisity  and  compos¬ 

ition  as  determined  by  a  mass  spectrometer  with  cryo  ion  source,  J. 

Geophys.  Res.  78:8296-8304. 


DURATION  OF  DISTURBANCE 
(DAYS) 


Figure  16.  Summary  of  Disturbed  Atmosphere  Frequency  or  Occurrence  (geo¬ 
magnetic  disturbances  only)  for  Solar  Maximum  Years 


by  Trinks  and  Fricke  (1978).^  In  Figure  17,  these  results  arc  plotted  as  ratios  to 
the  Ar  density.  At  night  it  is  found  that  the  C02/Ar  ratio  follows  the  expected 
diffusive  equilibrium  profile,  while  during  the  day  the  ratio  decreases  from  the 
diffusive  profile  by  a  factor  of  about  3  at  140  km.  This  result  is  in  agreement 
with  the  loss  of  CO2  by  solar  dissociation  and  reaction  with  O  ions.  The  fact 
that  the  CO2  density  responds  in  a  similar  way  to  the  Ar  density  makes  it  possi¬ 
ble  to  infer  the  behavior  of  fluctuations  in  the  CO2  density  from  measurements  of 
the  fluctuations  in  the  Ar  density,  especially  at  higher  altitudes. 

A  midlatitude  model  of  the  ozone  density  constructed  by  A.  J.  Kruger  and 
14 

R.A.  Minzner  was  presented  in  USSA  76.  Figure  18  (a  and  b)  shows  the  re¬ 
sultant  Og  concentration,  as  well  as  the  mass  mixing  ratio.  Expected  vari¬ 
ations  of  the  density  and  mixing  ratio  of  are  also  shown  in  these  figures.  From 
Figure  18  (a)  it  is  seen  that  at  high  altitude  the  concentration  varies  by  as  much 

as  a  factor  of  three.  During  the  Aladdin  program  in  1974,  measurements  of  the 

15 

Og  density  (Weeks  et  al,  1978  )  were  made  at  altitudes  from  19  to  108  km  using 

four  different  techniques.  Results  of  these  measurements  are  shown  in  Figure  19 
(a,  b,  and  c)  for  three  separate  altitude  regions.  These  results  are  in  basic 

14.  Kruger,  A.  J.  and  Minzner,  R.A,  (1976)  A  midlatitude  ozone  model  for  the 

1976  U.S.  standard  atmosphere,  J.  Geophys.  Res.  8 1:4477 -448 1. 

15.  Weeks.  L.H.,  Good,  R.E.,  Randhawa,  J.S.  and  Trinks,  H.  (1978)Ozone 

measurements  in  the  stratosphere,  mesosphere,  and  lower  thermosphere 
during  Aladdin  74,  J.  Geophys.  Res.  63:976-982. 
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Figure  17.  Calculated  and  Measured  n(C02/n(Ar)  Ratio 
for  Two  Background  Atmospheres  and  Three  Solar  Zenith 
Angles.  Data  are  from  Trinks  and  Fricke  (1978)®  (solid 
circles),  Offermann  and  Grossmann  (1973)^^(6quares) 
and  from  Philbrick  et  al,  (1971).  16  (Source  -  Trinks  and 
Fricke  (1978)9] 


agreement  with  the  model  presented  in  Figure  18  (a  and  b).  In  Figures  19  (a,  b, 
and  c)  these  results  are  also  compared  to  other  models  of  the  Og  density.  These 
measurements  reveal  a  structured  density  profile  at  altitudes  above  75  km,  par¬ 
ticularly  for  measurements  made  with  the  airglow  photometer.  These  variations 

16.  Philbrick,  C.R.,  Faucher,  G.A. ,  and  Wlodyka,  H.A.  (1971)  Neutral  Com¬ 
position  Measurements  of  the  Mesosphere  and  Lower  Thermosphere, 
A5’CrL-71-0602,  ad  739169.  - 


46 


ALTITUDE  (km) 


Figure  19a.  Chemiluminescent  Sonde  O3  Data  Com¬ 
pared  With  UV  Results,  ECC  Sonde,  and  Models 
(Weeks  et  al,  19781^) 
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Figure  19b.  Comparison  ol  UV  and  Airglow  Techniques  lor  O3  Uetei  mination 
(Weeks  et  al,  1978l'>) 
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Figure  19c.  Comparison  of  Airglow  and  Mass  Spectrometer  Re¬ 
sults  for  O3  Determination  With  Some  Typical  Models 
(Weeks  et  al,  19781^) 


appear  to  have  similar  scale  sizes  and  amplitudes  as  those  found  in  the  atomic 

oxygen  profiles  that  are  discussed  in  Chapter  9. 

Using  resonantly  scattered  light  from  thermospheric  NO,  measurements  made 

aboard  the  Atmosphere  Explorer  C  satellite  have  been  used  to  study  the  dependence 

of  the  NO  density  on  latitude,  longitude,  and  magnetic  activity.  These  results 

1 7 

have  been  analyzed  by  T.  E.  Cravens  and  A.  I.  Stewart  (1978)  to  determine  the 

NO  density  at  an  altitude  of  105  km.  They  find  that  near  the  equator  the  density 

7  -3 

at  the  peak  of  the  NO  concentration  is  about  2X10  cm  and  varies  little  with 
longitude  or  magnetic  activity,  except  during  major  geomagnetic  storms.  As 
shown  in  Figure  20,  the  peak  NO  densities  are  ;  pically  two  to  three  times  larger 

1  '7 

and  more  variable  at  higher  latitudes.  Cravens  and  Stewart  found  an  increase 
in  the  NO  density  of  as  much  as  30  percent  at  a  latitude  of  40°N  and  S.  In  addi¬ 
tion  to  these  effects,  they  found  that  the  NO  density  at  these  latitudes  doubles 
during  times  of  high  geomagnetic  activity  at  (K  (i).  They  also  studied  the 

17.  Cravens,  T.  E.  and  Stewart,  A,  I.  (1978)  Global  morphology  of  nitric  oxide  in 
the  lower  E  region,  J,  Geophys.  Res.  83:2446-2452. 
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diurnal  and  seasonal  effects  at  low  latitudes  and  found  that  NO  increases  at  the 
equator  by  a  factor  of  2  between  sunrise  and  early  afternoon  with  a  small  decline 
near  sunset  and  that  the  density  of  NO  is  higher  during  the  summer  than  during 
the  winter. 
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4.  Effects  of  Atmospheric  Variability 
on  Infored  Radiance  Prediction 


T.C.  Defsges 
Vuidyne,  Inc. 
Burlin^on,  Mags. 


The  infrared  radiance  model  described  in  Chapter  2  accepts  as  input  data,  molecu¬ 
lar  concentrations,  and  temperatures  at  fixed  altitudes.  The  computed  radiances 
depend  on  the  altitude  profiles  of  species  concentrations  and  temperature.  Vari¬ 
ability  in  the  predicted  radiances  results  from  the  variability  in  these  atmospheric 
properties  due  to  changing  geophysical  conditions.  In  this  chapter  we  present 
some  preliminary  results  that  show  the  effects  on  the  infrared  radiance  of  atmos¬ 
pheric  variations  as  described  in  Chapter  3. 

For  a  radiator  such  as  CO2,  which  is  chemically  stable  except  for  photodisso¬ 
ciation  by  far  ultraviolet  sunlight  at  high  altitudes  and  by  ion-molecule  reactions 
with  in  the  E  and  F  region,  temperature  is  the  controlling  factor  in  radiance 
calculations  below  100  km.  Above  100  km,  where  other  factors  influence  concen¬ 
trations,  it  is  still  the  radiation  received  from  lower  altitudes  and  the  sun  that 
controls  the  emitted  infrared.  Given  the  surface  pressure  and  concentration,  the 
atmospheric  temperature  profiles  determine  total  atmospheric  density  and  hence 
CO2  concentrations.  Up  to  30  or  40  km,  collisional  excitation  and  deexcitation  are 
oufficiently  rapid  to  maintain  vibrational  levels  of  COg  in  thermodynamic  equilibrium 
and  local  temperatures  determine  the  radiation  field.  Above  100  km,  collisional 
energy  exchange  is  much  slower  than  radiation  exchange.  The  radiation  from  the 
sun  and  lower  atmosphere  determines  the  radiation  absorbed  and  emitted  by  each 
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molecule.  Between  40  and  100  km,  radiation  and  collisions  compete  in  determin¬ 
ing  the  IR  emission.  Thus,  in  this  region,  the  entire  temperature  altitude  profile 
may  be  necessary  to  compute  radiation  from  molecules  at  a  given  altitude. 

Some  of  the  considerations  outlined  above  are  illustrated  in  Figure  1.  This 
figure  shows  total  band  radiance  in  the  CO2  01101-00001  (or  0110  -  0000)  15.  0  jjm 
vibration-rotation  band.  Four  different  nighttime  conditions  are  illustrated. 

Three  of  these  use  the  Spring/ Fall  Midlatitude  Model  from  the  USSA  7  6  for  alti¬ 
tudes  below  120  km.  The  600  K  exospheric  model,  which  represents  quiet  condi¬ 
tions  near  solar  activity  minimum,  is  further  assumed  to  have  a  COg  homopause 
(or  mictopause)  altitude  of  100  km.  The  2100  K  exosphere  models,  typical  of 
quiet  conditions  near  solar  activity  maximum  use  homopause  altitudes  of  110  and 
120  km.  The  fourth  set  uses  a  60°N  latitude  January  atmosphere  below  120  km, 
with  a  2100  K  exospheric  temperature  atmosphere  above  120  km.  The  homopause 
altitude  is  120  km.  Although  the  CO^  density  above  12  0  km  for  this  set  is  less  than 
that  of  the  corresponding  spring/fall  model,  radiance  values  are  higher  because 
the  upper  stratosphere  is  warmer  in  this  mode  and  the  "earthshine"  effect  (really 
"stratosphere  shine"  for  this  optically  thick  band)  is  stronger.  The  different 
radiance  values  for  the  spring/fall  models  are  due  almost  entirely  to  differing 
CO2  concentrations  at  high  altitude.  A  change  in  homopause  altitude  from  110  to 
120  km  approximately  doubles  concentrations  at  higher  altitudes,  while  varying 
exospheric  temperatures  produce  much  greater  variations  in  concentrations  and 
hence  radiation. 

The  USSA  76  supplements  provide  a  limited  set  of  reference  atmospheres  for 
illustrating  global  and  temporal  changes  in  atmospheric  conditions.  This  set  has 
the  advantage  that  consistent  atmospheres  are  given  from  the  ground  up  to  lOO  km. 
In  matching  the  high  altitude  parts  of  the  models  with  the  low  altitude  parts  there 
is  no  provision  for  possible  temperature  and  concentration  changes  below  120  km 
when  exospheric  temperatures  are  varied.  There  is  no  provision  for  including  the 
effects  of  magnetic  activity  (storms,  ring  current,  etc.  )  on  atmospheric  density, 
composition,  or  temperature.  As  discussed  in  Chapter  3  the  most  up-to-date 
model  atmospheres  that  incorporate  such  effects  are  the  J77  and  MSIS  models. 
These  have  the  disadvantage  that  they  are  limited  to  altitudes  above  PO  km  and 
cannot  be  expected  to  be  accurate  models  below  about  14  0  km.  the  lower  altitude 
limit  for  long-lived  satellites,  the  ultimate  source  of  the  data  used  in  constructing 
these  models. 

In  the  calculations  discussed  below,  some  of  the  problems  inherent  in  matching 
the  J77  Model  to  a  lower  altitude  model  have  been  avoided  by  selecting  models  from 
Cole  and  Kantor's  1978  Air  Force  Reference  Atmospheres  that  agree  closely  with 
the  J77  Model's  90  km  boundary  conditions.  Most  of  the  sample  radiances  have 
used  the  September  60°  N  latitude  model,  so  that  some  of  these  calculations  may 


Figure  1.  Vertical  (upward)  Viewing  of  Radiances  in  DllOl-oooDi  (.'O^ 
15  Band  for  Typical  Atmospheric  Models 


also  be  appropriate  to  analysis  of  the  AFGL  SPIRE  probe  data.  The  entire  J77 
model  has  not  been  used  but  the  temperature  was  taken  as  the  main  input  to  t'^e 
radiance  model  and  integration  is  performed  upward  from  the  .Air  Force  Reference 
Atmosphere  (Cole  and  Kantor,  1978^)  to  obtain  densities.  .At  9o  km,  the  total  mass 

density  from  the  60°  N  latitude  model  is  3.  33  X  lO  ^  kg-m  ^  as  compared  with 

*"6  *3 

3.43  X  10  kg-m  for  the  J77  density.  In  addition,  the  Jacchia  model  assumes  a 
homopause  at  100  km.  Here  the  effective  homopausc  level  for  carbon  dioxide  is 
taken  as  1 10  km. 

Eight  daytime  atmospheric  conditions  for  altitudes  above  90  km  have  been 
selected  for  use  with  the  September  60°  N  Cole  and  Kantor  model.  ^  These  condi¬ 
tions  are  summarized  in  Table  1. 

The  cases  with  solar  flux  300  (the  units  of  solar  flux  are  10  “  W  m  “  Hz  ) 

were  computed  at  7  5°  N  latitude  to  obtain  a  larger  effect  from  the  magnetic 


1.  Cole,  A.*  and  Kantor,  A.J.  (1978)  Air  Force  Reference  Atmospheres, 
AFGL-lK-78-0051,  AD  A058  505. 


T,  iC  1.  Atmospheric  Conditions  Used  for  Calculations  of  Infrared  Radiances 


Case 

Solar  Flux 

Exospheric  Temperature,  K 

A 

50 

0 

626.  5 

B 

100 

0 

880.  3 

C 

150 

0 

1081.2 

D 

2  00 

0 

1255.  0 

E 

250 

0 

1411.  5 

F 

300 

0 

1555.  7 

G 

300 

5 

1829.4 

H 

300 

8 

2176.  8 

disturbances.  The  60°  latitude  low  altitude  model  was  used  in  all  cases.  Plots 
of  band  radiances  presented  here  show  only  conditions  A,  C,  F,  G,  and  H. 

The  radiance  from  the  three  CO2  bands  with  band  center  at  15.  0  tim  are  added 
together  in  Figure  2.  The  01101-00001  transition  is  the  major  contribution  to  the 
15  Pm  bands  and  is  optically  thick  in  limb  viewing  up  to  a  tangent  height  of  about 
110  km.  For  this  transition  the  conditions  with  the  warmest  thermospheres  have 
sufficiently  more  molecules  at  high  altitudes  so  that  their  radiance  values  ai'e 
noticeably  higher  at  tangent  heights  as  low  as  60  km,  although  the  composition  is 
identical  for  all  cases  at  altitudes  below  90  km.  The  other  bands  do  not  have 
noticeable  differences  below  tangent  height  80  km.  Figures  3  and  4  show  total  band 
radiances  for  the  strongest  offcenter  vibrational-rotational  bands  included  in  the 
model  that  contribute  to  the  CO,  15  pm  band  complex.  Figures  5  through  7  show 
the  spectral  radiance  from  the  15  pm  CO2  bands  at  tangent  heights  100,  140,  and 
200  km,  with  2  percent  spectral  resolution.  At  100  km  the  cases  with  no  magnetic 
effects  have  almost  identical  radiance  curves,  while  the  cases  with  magnetic  effects 
G  and  H  are  brighter. 

Figure  8  shows  the  sum  of  three  bands  that  radiate  at  4.3  pm  for  cases  A,  C, 

F,  G,  and  H.  The  00011-00001  transition  band  is  optically  thick  in  limb  viewing 
to  about  12  0  km  tangent  height.  The  fluorescent  bands  are  brighter  than  the  funda¬ 
mental  between  60  and  100  km,  but  rapidly  become  less  important  at  higher  alti¬ 
tudes.  Figures  9  through  11  present  daytime  spectra  at  4.  3  pm  for  tangent  heights 
100,  140,  and  200  km.  As  was  the  case  at  15  pm,  the  cases  without  magnetic 
disturbances  are  very  close  together  at  100  km  tangent  height.  At  higher  altitudes 
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Figure  2.  Limb  Viewing  Radiance  Values  for  CO,  15  wm  Band 
Complex,  Cases  A,  C,  F,  G,  and  H 


density  differences  lead  to  separation  of  the  urves.  Figure  12  shows  the  radiance 
values  foi  the  10012-00001  CO2  band  at  2.77  /jin. 

Figures  13  and  14  present  radiances  from  two  important  bands  of  ^2^^'  ^l^®sc 
bands  ai'e  tire  Oi0-00(i  0.27  rrm  and  the  001-000  2.63  urn  bands.  Figures  1.5  and  lb 
compare  NO  band  radiances.  Figures  16  thr.jugh  18  give  5.4  ^m  NO  spectral 
radiances  at  tangent  heights  of  100,  140,  and  200  ian.  NO  concentrations  are  iden¬ 
tical  for  all  cases,  as  are  atomic  oxygen  concentrations.  The  brightness  differ¬ 
ences  below  140  km  tangent  height  are  due  to  tempr  rafurc  differences  that  control 
the  rate  of  excitation  by  collision  witli  atomic  oxygen.  The  O^  band  radiances  at 
0.  6  pm  are  shown  m  Figure  20. 

The  low  altitude  atmospheric  model.s  available  (see  Chapter  3,  Section  2,  3) 
are  models  of  average  conditions,  monthly,  at  a  small  number  of  latitudes.  There 
is  no  model  rea.iily  available  that  can  follow  the  daily  variations  in  tropospheric 
weather  and  stratospheric  circulation  patterns.  It  is  possible  to  simulate  vari¬ 
ability  by  adding  to  the  model's  smooth  temper  ature  profile  random  temperature 
excursions.  A  number  of  different  profiles  based  on  a  single  temperature  profile 
could  then  be  used  to  simulate  viewing  paths  through  a  varying  atmosphere. 
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LIMB  VIEWING  BANG  RflOIfiNCE.  WATTS/ ( CM»»2-STER ) 

Figure  3,  Limb  Viewing  Radiance  Values  for  the  10002-01101 
CO2  16.2  Mm  Band,  Cases  A,  C,  P,  G,  and  H 


The  matching  of  low  altitude  atmospheric  models  to  high  altitude  models  such 
as  MSIS  and  J77  requires  a  procedure  to  integrate  downward  from,  for  instance, 

150  km,  where  the  model  is  reasonably  accurate,  through  tlie  region  between  120 
and  90  km,  where  temperature  and  composition  are  not  known  well,  to  the  top  of 
the  low  altitude  model.  The  problem  is  in  determining  a  reasonable  temperature 
profile  to  bridge  the  gap  between  90  and  150  km.  In  each  case,  several  trials  will 
be  required  and  criteria  must  be  developed  to  make  the  downward  integration  fol¬ 
low  a  temperature  profile  that  is  in  accord  with  what  is  known  about  this  region. 

The  general  technique  is  available  in  the  algoriUims  used  in  analyzing  falling  sphere 
data. 
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Line  VIEWING  BAND  RADIANCE.  HATTS/t Cn»»2-STER ) 

Figure  4.  Limb  Viewing  Radiance  Values  for  the  10001-01101 
COg  13.  9  nm  Band,  Cases  A,  C,  F,  G,  and  H 
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Figure  7.  Limb  Viewing  Spectral  Radiances  Near  15 
Tangent  Height 
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Figure  9.  LimbViewing  Spectral  Radiances  Near4.3pin, 
100  km  Tangent  Height 
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Figure  10.  LimbVlewing  Spectral  Radiances  Near  4.3iiTn 
140  km  Tangent  Height 
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Figure  11.  LimbViewing  Spectral  Radiances  Near  4.3pm 
200  km  Tangent  Height 
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Figure  12.  Limb  Viewing  Radiance  Values  for  the  10012-00001 
CO2  2.77  »im  Band,  C.ases  A,  C,  F,  G,  and  H 
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Figure  13.  Limb  Viewing  Radiance  Values  for  the  010-100  H„0 
6.27  (tm  Band,  Cases  A,  C,  F,  G,  and  H 
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Figure  19.  Limb  Viewing  Radiances  for  NO  2-0  Overtone  Band, 
Cases  A,  C,  F,  G.  and  H 
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1.  INTRODUCTION 


In  addition  to  the  long-term  variations  of  the  atmospheric  species  concentra¬ 
tions  and  temperature  described  in  Chapter  3.  the  atmosphere  is  subject  to  vari¬ 
ations  that  have  shorter  periods  and  smaller  spatial  wavelengths.  These  variations 
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in  species  concentrations  and  temperature  can  result  in  fluctuations  of  the  result¬ 
ant  infrared  Intensity  that  can  effect  the  interpretation  of  infrared  measurements. 
In  this  chapter  we  review  information  concerning  these  variations  and  draw  some 
conclusions  as  to  the  magnitude  and  size  of  the  spatial  structures  as  functions  of 
altitude. 

A  preliminary  model  for  the  calculation  of  the  noise  (that  is,  the  clutter) 
created  in  infrared  measuring  systems  due  to  the  effects  of  random  spatial  struc¬ 
tures  on  the  radiant  intensity  is  also  discussed.  In  addition,  we  shall  briefly  dis¬ 
cuss  auroral  emissions  and  a  model  showing  how  spatial  fluctuations  in  TV  all-sky 
pl'-'l.'.netric  measurements  of  visible  aurora  radiation  can  be  used  to  generate 
spatial  power  spectral  densities  of  the  auroral  infrared  radiance. 


2.  SOURCES  OF  ATMOSPHERIC  STRUCTURE 

Dynamic  effects  such  as  tides  and  gravity  waves  can  perturb  a  "quiet  static" 
atmosphere.  The  sources  of  these  dynamic  effects  can  be  associated  with  a  very 
large  variety  of  phenomenon,  such  as  absorption  of  solar  radiation  by  atmospheric 
species,  solar  heating  of  the  surface,  sunspot  variations,  magnetospheric  sub¬ 
storms,  particle  events,  severe  tropospheric  disturbances,  planetary  waves,  wind- 
induced  waves  from  interaction  with  orographic  features,  earthquakes,  and  volcan¬ 
ic  eruptions  or  artificially  produced  phenomenon  such  as  bomb  explosions. 

2.1  Tidal  Effects 

Large-scale  tidal  oscillations  of  the  atmosphere  are  produced  mainly  by  the 
thermal  action  resulting  from  the  absorption  of  solar  radiation  in  the  stratosphere 
and  lower  mesosphere.  That  the  very  strong  semidiurnal  component  in  the  baro¬ 
metric  variation  was  due  mainly  to  the  sun's  heating  was  first  recognized  by 
Laplace.  The  present  inviscid  tidal  theory  is  based  on  the  classical  work  of 
Siebert  (1954)  and  Sen  and  White,  (1955),  who  were  the  first  to  consider  atmos¬ 
pheric  heating  rather  than  heating  by  upward  eddy  conduction  from  the  ground. 
Siebert  found  that  absorption  by  H2O  can  account  for  1/3  of  the  large  semidiurnal 
component  while  absorption  by  Og  can  account  for  the  remaining  two  thirds.  The 
large  semidiurnal  component  of  pressure  at  the  ground  level,  relative  to  the  di¬ 
urnal  component,  and  the  higher  values  of  the  diurnal  component  in  the  middle  and 
upper  atmosphere  are  explained  on  the  basis  that  heating  from  absorption  does 
not  reach  ground. 

1.  Siebert,  M,  (1954)  Zur  theoric  der  thermischen  erregung  gezeitenaitiger 

schwingunger  der  erdatmosphare.  Z.  Naturwissenschaften  41:446. 

2.  Sen,  H.  K.  and  White,  M.  L.  (1955)  Thermal  and  gravitational  excitation  of 

atmospheric  oscillations,  J.  Geophys.  Res.  60:483-485. 
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For  a  basic  treatment  of  tidal  wave  theory  the  reader  is  referred  to  Chapman 
and  Lindzen  (1970).  Recent  theoretical  advances  in  tidal  wave  theory  deal  mainly 
with  the  addition  of  such  effects  as  mean  zonal  winds,  meridional  temperature 
gradients,  composition  variations,  and  hydromagnetic  coupling,  as  well  as  dissi¬ 
pative  effects  such  as  molecular  viscosity,  thermal  conductivity,  radiative  damping, 
and  ion  drag.  For  a  review  that  considers  some  of  these  recent  advances  see 
Forbes  and  Garrett  (1978).  ^  The  empirical  classification  of  the  observed  tides, 
as  well  as  the  theoretical  treatment,  is  generally  done  in  the  form  of  harmonic 
analysis.  However,  owing  to  weather  and  seasonal  changes,  the  dependent  vari¬ 
ables,  such  as  pressure  or  temperature,  are  not  quite  cyclic  on  a  dav-to-daj  oasis. 
In  the  thermosphere  where  ion  drag,  viscosity,  and  heat  conduction  are  important, 
there  is  no  such  thing  as  a  pure  tidal  mode. 

The  tidal  effects  on  an  otherwise  unperturbed  atmosphere  produce  inhomogene¬ 
ities  in  optical  emission  and  composition  of  various  atmospheric  species,  some  of 
which  can  take  part  in  reactions  involving  infrared  radiation  bands.  The  effects 
on  the  emission  of  the  green  line  of  atomic  oxygen  (557.  7  nm)  were  considered  by 
Forbes  and  Geller  (1972)^  and  Petitdidier  and  Teitelbaum  (1.977).  ^  Assuming  the 
Chapman  three -body  recombination  mechanism  for  molecular  oxygen  and  chemical 
equilibrium,  the  green  line  intensity  is  proportional  to  the  cube  of  the  atomic  oxy- 

5 

gen  density.  Forbes  and  Geller  (1972)'  determined  the  variation  of  (O)  with  lunar 
tide  from  the  continuity  equation  neglecting  production  and  loss.  They  concluded 
that  the  lunar  tide  is  sufficient  to  produce  a  variation  in  the  green  line  airglow 
similar  to  that  observed.  Petitdidier  and  Teitelbaum  (1977)*^  considered  the  ef¬ 
fects  of  solar  thermal  tides.  They  assumed  that  the  composition  of  the  abnosphere 
■.  emains  unchanged  under  the  influence  of  the  tidal  wave  and  used  a  model  that 
assumes  a  sinusoidal  variation  of  the  relative  density  in  terms  of  time  and  space. 
Using  a  vertical  wavelength  of  45  km  and  a  period  of  12  h,  that  is,  a  semidiurnal 

tide,  Petitdidier  and  Teitelbaum  obtained  intensity  ratios  (I  .  /I  )  in  fair  agree- 

^  I'D  irr  iD 

ment  with  observed  data.  They  were  also  able  to  fit  observed  OH  experimental 
data  with  a  semidiurnal  tide.  This  analysis  did  not  give  a  very  sensitive  test  to 
the  tidal  wave  theory  since  they  compared  only  the  total  column  intensity  ratios. 

The  fc.rmation  of  OH  from  the  reaction 


3.  Chapman,  S. ,  and  Lindzen,  R.S.  (1970)  Atmospheric  Tides,  D.  Reidel, 

Hingham,  Mass. 

4,  Forbes,  J.  M.  and  Garrett,  H.  B.  (1978)  Solar  diurnal  tide  in  the  thermosphere, 

J.  Atmos.  Set.  35: 148-208. 

5,  Forbes,  J.  M.  and  Geller,  M.A.  (1972)  Lunar  semidiurnal  variation  in  OI 

(5577A)  nightglow,  J.  Geophys.  Res.  77:2942-2947. 

6.  Petitdidier,  M.  and  Teitelbaum,  H.  (1977)  Lower  thermosphere  emissions  and 

tides.  Planet  Space  Sci.  25:711-721, 
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(1) 


H  +  Og  —  OH  +  O2 

is  dependent  on  the  concentrations  of  the  minor  species  H  and  that  are  influ¬ 
enced  by  the  tidal  wave.  The  assumption  that  the  atmospheric  composition  is  not 
influenced  by  the  tidal  wave  is  only  true  for  major  species. 

Diurnal  variations  in  the  concentrations  of  H,  He,  O,  O2,  N2,  and  Ar  due  to 

energy  sources  in  the  lower  atmosphere  and,  in  particular  the  EUV,  have  been 

7  8 

considered  by  Harris  and  Mayr  (1975,  1977  ).  Their  calculations  are  based  on 

a  two-component  model  in  which  the  two  components  are  the  major  species  N2  and 
O.  This  treatment  assumes,  that  while  the  major  components  may  impart  momen¬ 
tum  to  the  minor  components,  the  minor  components  do  not  physically  or  chemi¬ 
cally  affect  the  major  components  in  any  way.  Their  results  show  that  the  most 
important  contribution  to  the  diurnal  variation  is  the  thermal  expansion  followed 
by  wind-induced  diffusion  that,  in  fact,  dominates  the  diurnal  variation  of  He  and 

prevails  strongly  in  the  variations  of  O,  0„,  and  Ar  below  200  km. 

Q  *  4 

Forbes  (1978)  and  Forbes  and  Garrett  (1978)  used  a  single  component  Ng  as 
the  major  species  and  all  other  constituents,  that  include  O,  O2,  Ar,  He,  and  H, 
as  minor  species.  The  minor  species  were  coupled  to  the  major  species  through 
collision.  Special  emphasis  is  placed  in  their  calculations  on  the  diurnal  and 
semidiurnal  variations  at  sunspot  maximum  and  minimum.  Their  results  show 
that,  between  the  equator  and  the  midlatitude,  the  diurnal  tidal  variation  in  the 
thermospheric  composition  depends,  at  a  given  height,  more  strongly  on  the  solar 
cycle  than  on  the  latitudinal  position.  The  opposite  is  true  for  the  semidiurnal 
tidal  variations,  where  they  find  a  major  dependence  with  latitude  in  the  semidi¬ 
urnal  amplitude  and  phase.  These  effects  are  presently  not  modeled  in  the  J77 
and  MSIS  modeis.  Differences  in  the  total  density  variations  between  equatorial 
and  middle  latitudes  are  predicted  by  this  theory  (Forbes  and  Marcos,  197  9)^*'  to 
be  of  the  order  of  20  -  50  percent  in  amplitude  and  3  -  6  h  in  phase  for  the  semi¬ 
diurnal  component,  0-20  percent  in  amplitude  and  0-6  h  in  phase  for  the  diurnal 
component,  and  15  km  in  the  altitude  of  transition  from  semidiurnal  to  diurnal 


7.  Harris,  I,  and  Mayr,  H.G.  (1975)  Diurnal  variations  in  the  thermosphere, 

J.  Geophys.  Res.  80:3925-3933. 

8.  Mayr,  H.G.  and  Harris,  I.  (1977)  Diurnal  tide,  composition,  and  winds, 
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predominance.  Experimental  data  that  supports  their  results  is  presented  in  the 
next  section. 


2.2  Atinospherir  Gravity  Waves 

Gravity  waves  are  atmospheric  waves  with  frequencies  lower  than  the  natural 
oscillation  of  the  atmosphere,  the  Brunt -Vaisala  frequency,  that  has  a  period 
ranging  from  less  than  5  min  to  15  min.  For  a  gravity  wave  the  pressure  gradi¬ 
ent,  which  may  be  along  any  direction,  is  modified  by  the  influence  of  gravitational 
and  buoyant  forces  that  are  along  the  vertical  direction.  Since  the  frequency  of  the 
wave  is  less  than  the  Brunt-Vaisala  frequency,  the  influence  of  the  gravity  and 
buoyancy  forces  begin  to  manifest  themselves  and  the  wave  is  no  longer  in  the 
direction  of  the  pressure  gradient.  For  a  free  gravity  wave  in  an  isothermal 
atmosphere  an  upward  transfer  of  energy  (group  velocity)  is  always  accompanied 
by  a  downward  phase  propagation  and  vice  versa  (Hines,  1900). 

As  far  as  introducing  variabilities  and  inhomogeneties,  there  are  three  main 
features  that  distinguish  gravity  waves  from  the  tides.  These  features  are: 

1)  gravity  waves  can  often  have  length  scales  of  the  order  of  0  to  8  km  (for  ex¬ 
ample,  vertical  wavelength)  that  are  considerably  less  than  the  length  scales  of 
tidal  waves,  2)  gravity  waves  can  have  time  scales  of  the  order  of  10  min  to  two 
or  three  hours  that  are  considerably  less  than  the  time  scales  of  dirunal  or  semi¬ 
diurnal  tides,  and  3)  gravity  waves  occur  randomly,  unlike  the  tides  that  occur 
regularly  due  to  solar  heating  of  the  atmosphere. 

The  effects  of  gravity  waves  were  initially  detected  in  mici’obaragraphs  at  the 
ground  level,  airglow  measurements  from  various  altitudes  in  the  atmosphere, 
and  ionospheric  radar  measurements  from  higher  altitudes.  Recent  techniques 
allow  higli  resolution  in-situ  measurements  of  the  effects  of  gravity  waves.  The 

earliest  observations  were  made  indirectly  through  their  effects  on  the  ionosphere 
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called  traveling  ionosphere  disturbances  (TID).  Georges  (1988)  characterized 
the  spectrum  of  observed  TID  into  two  distinct  categories  of  waves:  1)  the  large- 
scale  TID's  have  horizontal  phase  velocities  greater  than  33()-4()0  m/sec  and  hori¬ 
zontal  wavelengths  exceeding  1000  km  and  2)  the  medium-scale  TID's  have  hori¬ 
zontal  phase  velocities  commonly  between  100  and  250  m/sec  and  horizontal  wave¬ 
lengths  of  the  order  of  several  hundred  km's.  They  occur  much  more  frequently 

than  the  large-scale  TID's.  The  Significant  contribution  in  establishing  that  these 
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large-scale  TID's  are  caused  by  gravity  waves  was  made  by  Thome  (1988),  who 
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used  the  measurements  of  incoherent  radar  backscatter  from  different  altitudes 
in  the  ionosphere  as  a  function  of  time. 

Gravity  waves  have  also  been  detected  from  their  effects  on  airglow  emission. 
Many  of  the  airglow  emissions  (for  example,  O  557.7  nm  line,  OH,  etc.  )  have 
their  peak  intensities  well  below  the  peak  of  the  ionospheric  F-layer.  This  allows 
studies  to  be  made  over  larger  altitude  ranges  whenever  the  airglow  data  can  be 

combined  with  the  ionospheric  data  (for  example,  Dachs  (1968)^^  and  Armstrong 
15 

(1975)  ].  Airglow  data  showing  wave-like  structures  have  been  published  by 

numerous  people  including  Silverman  (1962),  Okuda  (1962),  and  Dachs  (1968). 
Dach's  data  (1968)  include  observations  on  the  O  lines  at  557.7  nm  and  630.0  nm 
as  well  as  simultaneous  observations  of  various  ionospheric  parameters  such  as 
the  altitude  of  the  F2  peak  h'F2  and  the  lowest  penetration  frequency  at  the  peak 
f®F2. 

The  basic  paper  on  gravity-wave  theory  is  due  to  Hines  (1960),  who  assumed 
an  inviscid  uniform  isothermal  atmosphere.  Even  with  such  drastic  assumptions 
this  treatment  provides  a  reasonably  good  first  approximation  for  gravity  waves 
with  low  horizontal  phase  velocities  (less  than  150  m/sec).  For  higher  horizontal 
phase  velocities  the  Ueatment  breaks  down  as  atmospheric  structure  begins  to 
have  an  effect.  Investigations  of  the  acoustic  gravity  equations  for  a  more  real¬ 
istic  atmosphere  show  the  presence  of  fully  guided  modes  [  Press  and  Harkrider 

(1962),  Pfeffer  and  Zarichny  (1963),  Harkrider  and  Wells  (1968)^'^)  as  well 
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as  partially  guided  or  leakage  modes  [Friedman  (1966),  Reddy  (1969)"  ].  The 
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inclusion  of  dissipation  can  produce  additional  partially  guided  modes  [  Tuan 
(1976),  Richmond  (1978)  ],  In  general,  inhomogeneities  can  scatter  gravity 

waves  and  hence,  may  (although  not  necessarily)  produce  additional  guided  modes. 

Both  for  the  case  of  the  tidal  waves  and  gravity  waves  the  theories  have  been 
mainly  tested  against  the  observed  temporal  variations  in  the  optical  emission. 

The  effect  of  gravity  waves  on  the  atmospheric  emission  of  the  O  630  nm  and  the 

2  5 

N  520  nm  lines  was  treated  theoretically  by  Porter.  Silverman,  and  Tuan  (1974) 
with  satisfactory  agreement  between  theory  and  experimental  data.  Spatial  inhomo¬ 
geneities  have  been  obseirved  in  OH  emission  by  Peterson  and  Kieffaber  [  Peterson 

P7 

and  Kieffaber  (1973),  Kieffaber  and  Peterson  (1978)  ].  They  found  well  de- 
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fined  OH  structures  to  occur  near  the  time  of  high  lunar  tide.  Peterson  (1979) 
suggested  that  these  structures  may  have  been  caused  by  gravity  wave  interactions 
in  the  OH  layer.  Gravity  wave  effects  on  any  optical  emission  depend  on  the 
atmospheric  minor  constituents  redistribution  caused  by  the  wave.  The  major 
constituents  are  assumed  to  move  with  the  wave  while  the  effects  of  the  wave  on 
the  minor  components  are  through  collisions  with  the  major  constituents. 

Chiu  and  Ching  (1978)  have  considered  the  general  features  that  a  Hines 
(1960)^^  gravity  wave  model  will  have  on  a  minor  atmospheric  component  whose 
density  distribution  is  in  the  form  of  a  single  layer.  Their  summary  of  the  conclu¬ 
sions  of  various  authors  are  1)  the  response  of  the  layer  to  the  gravity  wave  is 
sharper  at  the  bottom  side  than  at  the  top  side,  2)  there  is  a  180°  phase  change  in 
the  response  of  the  minor  constituent  in  going  from  the  bottom  to  the  top  side,  and 
3)  the  magnitude  of  the  response  depends  on  the  gradient  of  the  unperturbed  layer. 
The  sharper  the  gradient,  the  sharper  becomes  the  response. 
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Recently,  Hatfield  et  al  (1981)  studied  the  variability  of  the  OH  emission 
produced  by  gravity  waves.  Their  results  indicate  that  structure  in  the  emission 
profile  could  occur  if  the  vertical  wavelength  of  the  gravity  wave  is  less  than  the 
"width"  of  the  emission  profile.  This  corresponds  to  gravity  waves  with  horizon¬ 
tal  phase  velocities  less  than  33  m/sec.  Guided  modes  usually  have  much  greater 
horizontal  phase  velocities  and  hence,  much  larger  vertical  wavelength  and  can 
only  affect  the  integrated  columnar  emission  as  a  function  of  time.  Thus,  gravity 
waves  can  always  produce  temporal  inhomogeneities,  but  only  those  with  small 
horizontal  phase  velocities  and  vertical  wavelength  of  less  than  10  km  can  produce 
spatial  as  well  as  temporal  inhomogeneities. 

Further  studies  on  the  fluctuation  of  OH  and  0~(^E)  airglow  have  been  made 
31  1 

by  Weinstock  (1978).  He  found  that  for  0^(  L)  the  temperature  dependence  of 
the  reaction  rate  constant  has  a  major  influence  on  the  temporal  inhomogeneities. 

Spatial  and  temporal  inhomogeneities  can  be  generated  by  atmospheric  grav¬ 
ity  waves  if  the  energy  of  the  waves  is  dissipated  in  forming  a  layer  of  turbulence 
32 

[Hodges  (1967)  ].  In  some  cases  it  is  possible  to  show  that  large  amplitude 

gravity  waves  with  a  Richardson's  number  less  than  one-half  can  parametricallv 
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excite  Brunt-Vaisala  oscillations  (Tuan,  et  al,  1979  ).  This  means  that  temporal 

inhomogeneities  of  the  relatively  short  time  scales  of  Brunt-Vaisala  oscillations 
(of  the  order  of  4-12  min)  are  superimposed  on  inhomogeneities  of  the  larger  time 
scales  (or  the  order  of  hours)  of  the  gravity  wave  oscillations. 

Gravity  waves  that  are  usually  present  in  the  middle  and  upper  atmosphere 
are  frequently  associated  with  sources  in  the  lower  atmosphere.  At  times  tliis 
wave  propagation  can  be  effectively  cut  off  from  the  upper  atmosphere  by  reflet'- 
tion  or  dissipation  processes  in  the  mesosphere.  However,  sources  in  the  upper 
atmosphere  also  exists.  An  indication  of  the  horizontal  and  vertical  wavelengths 
of  gravity  waves  in  the  mesosphere  and  the  lower  thermosphere  is  given  in  Fig¬ 
ure  1  (Philbrick  et  al,  1978  ).  Also  shown  are  the  viscous  dissipation  limits  for 
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Figure  1.  Wavelengths  of  Pi'opagating  Modes 
Expected  in  the  Mesosphere  and  Lower  Thermo¬ 
sphere  From  a  Reexamination  of  the  Analysis 
of  Hines^^  With  Inclusion  of  Eddy  Viscosity 
Contribution.  The  periods  in  min  are  shown  as 
constant  period  contours  (solid  lines).  The 
viscous  dissipation  limit  for  small-scale  waves 
at  60  and  110  km  are  shown  (dash  lines)  and  re¬ 
present  the  limits  of  the  permitted  spectrum 
from  effects  of  viscous  damping  (modes  lying 
to  the  left  and  below  these  curves  are  excluded). 
Modes  subject  to  reflection  back  toward  the 
ground  at  heights  of  54  and  79  km  are  indicated 
(dot  lines)  and  inodes  lying  to  the  right  of  these 
curves  cannot  proceed  from  ^e  lower  atmos¬ 
phere  (Philbrick  et  al,  1978’'^). 
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small-scale  waves  and  the  limits  for  modes  subject  to  reflection  back  toward  the 
ground.  From  Figure  1  we  see  that  gravity  waves  propagating  in  the  region  Irom 
fiO  to  110  km  altitude  will  have  wavelengths  larger  than  1  km  at  (lU  km  and  larger 
than  3  km  at  110  km. 

Gravity  waves  may  be  produced  by  a  number  of  mechanisms  including  orogra¬ 
phic  features  such  as  mountain  ranges,  tropospheric  disturbances  such  as  thunder¬ 
storms,  and  localized  pressure  gradients  produced  by  auroral  heating  in  the  meso¬ 
sphere  and  lower  thermosphere.  Studies  have  been  made  by  Davis  and  Da  Rosa 
(iPtiO)  on  the  frequency  of  occurrence  of  gravity  waves  produced  by  magnetic 
stui  iiis  .  The  observations  for  the  year  1907  showed  a  close  correlation  between 
the  occurrence  of  TID's  and  severe  magnetic  storms  I'or  smaller  storms 

there  was  far  less  correlation.  These  results  are  in  good  agreement  with  those 
12  13 

of  Georges  (1908)  and  Thome  (1968)  who  suggested  a  relationship  between 
enhanced  magnetic  activity  and  the  occurrence  of  large-scale  TID’s. 

In  general,  there  has  been  little  investigation  on  the  frequency  of  occurrence 
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of  gravity  waves  due  to  other  sources.  If  the  theory  as  advanced  by  Hodges  (1967) 
that  gravity  waves  are  an  important  factor  in  producing  upper  atmospheric  turbu¬ 
lence  is  correct,  then  the  frequency  of  occurrence  must  be  much  greater  than  that 
from  magnetic  storms. 

3.  .ATMOSPHERE  STRliCTURE  MEASUREMENTS 

3.1  Tidal  EfferU 

The  effects  produced  by  the  thermally  driven  tides  in  the  upper  atmosphere 

present  another  important  consideration.  Recently  Forbes  and  Marcos  (197  9)'*' 

analyzed  total  mass  density  data  from  the  miniature  electrostatic  accelerometer 

(Mesa)  experiment  on  the  AE-E  satellite  to  determine  diurnal  and  semidiurnal 

tidal  variations  in  the  lower  thermosphere  (150-245  km)  at  latitudes  less  than  20°. 

The  results  of  the  analysis  showed  that  structure  changes  from  predominantly 

semidiurnal  below  180  ±5  km  to  predominantly  diurnal  above  this  altitude,  with 

no  strong  dependence  on  season.  These  results  are  in  agreement  with  recent 

8  9 

theoretical  (Mayr  and  Harris,  1977,  Forbes,  1978  )  and  observational  (Newton 
37  38 

et  al,  1975  Sharp  et  al,  1978  )  investigations  (see  Figure  2).  The  analysis  by 

Forbes  and  Marcos  showed  a  substantial  summer/winter  difference  in  the  phase 
structure  of  the  semidiurnal  component  of  the  total  mass  density  variation.'*'  The 
results  indicate  that  the  combination  of  data  from  all  seasons  to  determine  an 


References  36  to  38  will  not  be  listed  here.  See  References,  page  114. 
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LOCAL  SOLAR  TIME  (hr) 


Figure  2.  Average  Ratio  of  the  Measured 
Density  to  the  Jacchia  (1971)  Mode)  Predic¬ 
tion  as  a  Function  of  Local  Solar  Time.  The 
local  time  parameter  for  tlie  model  is  taken 
to  be  local  midnight  in  all  cases,  but  all  other 
parameters  are  handled  in  the  manner  pre¬ 
scribed  in  the  model  formulation.  The  data 
are  binned  at  two-h  intervals.  Each  curve 
includes  all  altitude  increments  withii.  the 
2  0  km  band  centered  on  the  designated  alti¬ 
tude  (Sharp  et  al,  1978^^) 


annual  mean  semidiurnal  amplitude  can  lead  to  unrealistically  small  semidiurnal 
amplitudes,  due  to  phase  cancellation  effects. 

3.2  \lim>9|ilipric  (iravity  Waves 

A  few  mea.surements  of  the  total  atmospheric  density  liave  been  made  using 
rocketborne  falling  spheres  with  accelerometer.':  ll)at  are  capable  of  making  meas¬ 
urements  with  altitude  re.solution  of  ILiO  iii.  These  measurements  reveal  a  high 
degree  of  spatial  structure  in  the  measured  density  (sec  Figures  3  and  4).  Fig¬ 
ure  3  shows  an  example  of  the  wave  structure  observed  in  the  mesosphere  and 

lower  thermosphere  using  a  recently  developed  accelerometer  technique  (Philbrick 
34 

et  al,  1978  ).  Small-scale  wave  structures  that  ar  e  consistent  w  ith  gravity  waves 

(see  Figure  1)  are  observed  to  be  superimposed  on  longer  wavelength  oscillations. 
The  large-scale  waves  above  70  km  altitude  in  Figures  3  and  4  are  presumably  due 
to  a  combination  of  large-scale  gravity  waves  of  tidal  oscillations  or  possibly  arti¬ 
facts  due  to  the  atmospheric  model  used  for  comparison.  The  appearance  of  the 
small-scale  waves,  as  shown  in  Uiese  figures,  has  been  obsci'ved  on  all  of  the  half 
dozen  flights  that  could  have  detected  such  waves.  Thus,  they  are  expected  to  be 
present  under  most  conditions.  Figure  4b  shows  the  temperature  derived  on  the 
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Figure  3.  Density  Measurements  Made  on  2fi  February  1970  are  Shown  as  a 
Ratio  to  the  USSA  76  Standard  Model.  Short  wavelength  (2-4  km)  structure 
features  are  superimposed  on  longer-scale  waves.  The  scatter  of  the  meas¬ 
urement  points  at  higher  altitudes  is  due  to  the  telemetry  resolution  of  the 
data 

basis  of  the  density  profile  of  Figure  4a  and  shows  a  high  degree  of  vertical  struc¬ 
ture. 

Figure  5  indicates  the  temporal  dependence  of  the  vertical  structure  (Philbrick 
30 

et  al,  1980  '  ).  This  l.^jOire  shows  the  results  of  density  measurements  made  at 
Kwajalein  Atoll  for  two  sets  of  falling  sphere  flights.  The  larger  waves  are  seen 
to  be  nearly  stationary  in  space  (that  is,  not  propagating  vertically)  over  periods 
of  several  hours.  These  rncasurement.s  were  accompanied  by  simultaneous  meas¬ 
urements  of  the  wind  velocity  that  permits  a  cah'ulation  of  the  Richardson  number 


39.  Philbrick,  (  .(!.,  Murphy,  1  .  \..  /imii  i-rman.  S.  P.  ,  I'Utcher,  R..I.,  3r.  . 
and  Olsen,  H.O.  (1980)  Me.sospheric  ilensitv  variability,  Sp;o  c  He.scan  h 
XX:79-82.  —  - - 
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Figure  4a.  Hypersonic  Sphere  Density  Results  from  the  Lincoln  Laboratory, 
Analysis  of  the  ALCOR  Radar  Data  Compared  to  Model.  Note  the  scale  sizes  of 
the  vertical  structure  compared  to  the  scales  expected  based  on  arguments  pre¬ 
sented  in  Figure  1  (Philbrick  et  al,  1978^'*) 


R^.  The  results  of  this  calculation  are  also  shown  in  Figure  5.  These  plots  show 
the  scale  sizes  and  the  dynamically  unstable  regions  (RjS  0.25)  due  to  the  presence 
of  wind  shear.  Negative  values  of  the  Richardson  number  exists  in  regions  that  are 
statically  unstable  because  of  the  local  density  gradient. 

In  addition  to  these  total  density  measurements,  several  rocketborne  mass 
spectrometers  have  made  measurements  in  the  altitude  range  from  80  to  120  km. 
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Figure  4b.  Temperature  from  Lincoln  Laboratory,  Analysis  of  the  Hypersonic 
Sphere  Data  Compared  to  Model  (Philbrick  et  al,  1978^'*} 

34 

Figure  6  shows  an  example  of  these  measurements  (Philbrick  et  al,  inTS  ), 
where  vertical  wave  structures  are  apparent  in  the  species  profiles.  The  atomic 
oxygen  profile  reveals  a  substantial  amount  of  spatial  structure.  Similar  struc¬ 
ture  has  been  observed  in  most  flights.  Measurements  of  the  vertical  structure 
of  atomic  oxygen  made  at  the  White  Sands  Missile  Range  (WSMR)  using  resonant 
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Figure  5,  Density  Measurements  and  Richardson  Numbers  for  Falling  Sphere 
Flights  at  Kwajalein  Atoll  on  18  May  1977  and  5  April  1978.  (a)  The  density  meas- 

urements  are  shown  as  a  ratio  to  the  Cole  and  Kantor  tropical  atmosphere  for  four 
Robin  sphere  flights  and  one  accelerometer  instrumented  sphere  (AC-6),  (b)  The 

Richardson  numbers  calculated  from  the  temperature  and  wind  measurements  are 
shown  for  the  range  of  values  0  <  R,  <  1  with  the  negative  values  plotted  at  R.  =  0. 
(PhUbrick  et  al.  198039)  i  i 


Figure  6.  Species  Densities  vs  Altitude  Measured  by  the 
Mass  Spectrometer  on  23  January  1976  (Philbrick  et  al, 
1978-^‘*) 


scattering  of  light  indicates  similar  spatial  structure  (Figure  7,  Hewlett  et  al, 
1980^°), 

Results  obtained  from  a  dynamic  model  of  the  mesosphere  and  lower  thermo- 

41 

sphere  (Keneshea  et  al,  1979  )  are  compared  in  Figure  8  with  profile  measure¬ 

ments  obtained  on  20  November  1970,  This  model  requires  data  on  the  atmospher 
ic  composition  and  temperature  as  a  function  of  altitude,  as  well  as  the  molecular 
and  turbulent  diffusion  coefficients.  The  latter  are  determined  from  wind  trail 
measurements.  The  substantial  agreement  between  theory  and  measurements  in¬ 
dicates  the  importance  of  considering  the  atmospheric  dynamics.  Recent  improve 
ments  in  the  Schumann-Runge  flux  factor  would  produce  less  atomic  oxygen  in  the 
model  calculation  and  improve  the  agreement  between  the  model  and  the  measured 
density  of  this  species.  The  agreement  supports  the  conclusion  that  atmospheric 
dynamics  can  strongly  influence  the  density  profile  of  a  minor  species  such  as 
atomic  oxygen  in  the  80  to  100  km  region. 

40.  Howlett,  L,  C. ,  Baker,  K.D.  ,  Megill,  L.R..  Shaw,  A.  W.  and  Pendleton, 

W.R.  (1980)  Measurement  of  a  structured  profile  of  atomic  oxygen  in  the 
mesosphere  and  lower  thermosphere,  J.  Geophys.  Res.  85:1291-1298. 

41.  Keneshea,  T.  J.  ,  Zimmerman.  S.  P.  and  Philbrick,  C.  R.  (1979)  A  dynamic 

model  of  the  mesosphere  and  lower  thermosphere.  Planet.  Space  Sci. 

27:  385-401.  - 
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ATOMIC  OXYGEN  DENSITY  ( Atom*  cm"®) 

Figure  7a.  Atomic  Oxygen  Density  on 
Ascent  Showing  an  Estimate  of  Hel«^t)ve 
and  Absolute  Errors  (Three-sec  lun- 
ning  average  with  one  point  per  iec 
shown)  (Hewlett  et  al,  1980^0) 


COUNT  RATE  (s*e'') 


Figure  7b.  High -resolution  Ascent  Data 
Showing  Con.siderable  Structure  (One -sec 
running  average  with  one  point  per 
0.  1  sec  shown)  (Howlett  et  al,  1980^^ 
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Figure  8.  Comparison  of  Computed  Species  Profiles  for  Turbulent  Dif- 
fusivities  K  and  3K  (solid  dashed  cujrves)  With  Measured  Profiles  for  O 
(squares),  O2  (triangle),  N2  (circles),  Ar  (crosses),  and  02(^Ag)(dots) 


At  higher  altitudes,  satellite  measurements  show  wave -like  variations  that  are 
thought  to  be  due  to  gravity  waves.  In  Figure  9,  measurements  made  during  flights 
of  the  S3-4  satellite  (Mclssac,  1980'^^)  reveal  the  presence  of  these  waves  with 
horizontal  wavelengths  of  hundreds  of  km.  These  waves  are  frequently  observed 
at  high  altitudes  and  are  always  present  for  disturbed  periods  (Kp  >4).  The  struc¬ 
ture  is  strongly  correlated  with  the  auroral  region  and  probably  associated  with 
pressure  gradients  created  by  momentum  sources  or  localized  heating  in  the  lower 
thermosphere  or  upper  mesosphere. 

The  nature  of  these  waves  is  shown  more  clearly  in  the  measurements  made 

4  3 

using  the  Nace  instruments  on  flights  of  the  AE-C  satellite  (Reber  et  al,  1975  ). 

In  Figure  10  the  wave  structure  is  clearly  revealed.  The  authors  of  this  study 
43 

(Reber  et  al,  1975  )  concluded  from  their  analyses  of  the  spatial  structure  for  two 

orbits  that;  "This  structure  exhibits  (1)  scale  sizes  from  a  few  tens  of  kilometers 
to  several  hundreds  of  km,  (2)  molecular  nitrogen  density  amplitudes  of  the  order 
of  5-15%,  argon  amplitudes  about  half  of  Ng  amplitudes  and  generally  180°  out- 
of -phase  with  the  and  Ar,  .  .  The  neutral  data  are  qualitatively  consistent 

42.  Mclssac,  J.  P.  (1980)  Private  communication. 

43.  Reber,  C.  A.,  Hedin,  A.E.  ,  Pelz,  D.  T. ,  Potter,  W.  E..  and  Brace,  L.  H. 

(1975)  Phase  and  amplitude  relationships  of  wave  structure  observed  in  the 
lower  thermosphere,  J.  Geophys.  Res.  80:4570-4580. 


90 


AMPERES  AMPERES  AMPERES 


S3-4  NEUTRAL  DENSITY  MEASUREMENTS 
Particle  Flux  Accumulator  Experiment 
NEUTRAL  DENSITY  PHENOMENA  -  Atmospheric  Waves 
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Figure  9.  S3-4  Neutral  Density  Measurements  —  Particle  Flux  Accum¬ 
ulator  Experiment  Neutral  Density  Phenomena  —  Atmospheric  Waves 
(Mclssac, 


An 

n 


8  FEB  1974 


05  32  UT 


LONGITUDE  109<> 


LAI  6J*  M*  ‘.(1“  42”  14*  24'  22'  i /' 

alt  iOO  2bO  :on  l/b  |*.8  1 /b  2UO  2t.L'  .bio  km 


Figure  10.  The  Structural  Variations  of  Molecular  Nitrogen.  Argon,  Total  Ion 
Density,  and  Electron  Temperature.  The  relative  scale  sizes  along  the  satellite 
trajectory  are  shown,  and  again  the  vertical  lines  emphasize  the  phasing  between 
parameters.  The  entire  data  set  was  obtained  in  14.7  min  at  a  satellite  speed  of 
8.  6  km*^  (Reber  et  al,  197543) 


with  a  simple  picture  of  wave-like  density  enhancements  in  the  major  species 
nearly  in  phase  with  the  vertical  component  of  the  perturbation  velocity,  ..." 

Later  analysis  of  the  OSS  data  obtained  on  the  AE-C  satellite  (Potter  et  al, 

44 

1976  )  has  shown  the  amplitudes  of  the  waves  increase  at  the  high-latitude  auroral 

regions,  indicating  that  the  waves  were  apparently  generated  by  magnetic  storm 
energy  deposition  and  propagate  both  northward  and  southward  from  the  region  of 
the  auroral  electrojet.  An  example  of  the  OSS  results  are  shown  in  Figure  11.  It 
is  expected  that  similar  wave  structure  would  appear  in  the  infrared  radiation  due 
to  CO2  from  these  regions,  as  the  spatial  variations  in  the  102  that  result  from 
mass  motions  are  expected  to  be  similar  tc  'hat  of  Ar. 

The  momentum  source  is  one  of  the  features  of  the  thermosphere  that  mani¬ 
fests  itself  in  significant  high-latitude  structure  during  periods  of  enhanced  geo¬ 
magnetic  activity.  The  momentum  source  manifests  itself  by  a  pressure  gradient 
due  to  momentum  transfer  to  neutrals  from  ions  accelerated  by  electric  fields. 

44.  Potter,  W.E.,  Kaysei',  D.  C.  and  Mnuersberger.  K.  (1976)  Direct  measure¬ 
ments  of  neutral  wave  characteristics  in  the  thermosphere,  .1.  Geophys. 
Res.  81:  ,5002-.5ni2.  - 
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Figure  11.  The  Center  and  Bottom  Panels  Show  OSS  Neutral 
Number  Densities  for  Two  Elliptical  Orbits.  The  mass  de¬ 
pendent  wave  structure  is  evident,  as  is  the  out-of-phase  be¬ 
havior  of  He.  The  top  panel  shows  the  total  Mesa  mass  density 
from  orbit  629,  that  closely  follows  the  OSS  N2  number  density 
in  the  center  panel.  The  circled  pluses  shown  in  the  top  panel 
are  mass  densities  calculated  from  OSS  constituent  number 
densities.  Altitudes  are  given  in  km.  The  (0)  densities  re¬ 
present  the  total  atmosphere  oxygen.  The  universal  time  (UT) 
is  given  in  sec  after  the  start  of  the  day 
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Perturbations  in  the  structure  of  the  neutral  atmosphere  in  narrow  latitude  bands 

near  the  dayside  polar  cusp  have  been  measured  aboard  the  S3-1  and  Esro  4  satel- 

45 

lites  (Trinks  et  al,  1978a  ).  Examples  of  these  perturbations  are  shown  in  Fig¬ 

ures  12  to  14.  Figures  12  and  13  present  measurements  made  aboard  the  S3-1 
satellite.  These  measurements  show  fluctuation  amplitudes  of  as  much  as  a  fac¬ 
tor  of  two  in  species  concentrations  and  with  spatial  dimensions  of  approximately 
300  km.  Figure  14  presents  the  measurements  made  aboard  the  Esro  4  satellite 
and  indicate  similar  magnitudes  and  extents.  Trinks  et  al,  (lOVSb)"^^  discuss 
these  measurements  and  suggest  that  thermospheric  winds  driven  by  "momentum 
sources"  associated  with  ion  convection  are  the  predominant  cause  for  these  dis¬ 
turbances.  The  narrow  spatial  structure  and  in  phase  response  of  all  species  are 
characteristics  associated  with  the  momentum  source.  These  effects  are  generally 
superimposed  on  the  thermosphere  expansion  changes  caused  by  Joule  heating. 
Momentum  source  disturbances  probably  generate  waves  which  are  propagated  to 
lower  latitudes. 

3.3  lonoepherir  Response  to  Geomagnetic  Activity 

Ion  densities  have  been  measured  on  a  number  of  satellite  flights.  The  effects 
of  increased  geomagnetic  activity  on  the  ionic  concentration  can  be  seen  by  compar¬ 
ing  Figures  15  and  16.  Figure  15  shows  ion  species  concentrations  during  a  quiet 
period  and  displays  the  typical  structure  of  the  F  layers.  In  Figure  16  the  normal 
F2  layer  at  high  latitudes  is  severely  changed  during  a  period  of  high  geomagnetic 

activity  (K  =  6).  In  the  disturbed  region  there  is  a  large  amount  of  spatial  struc- 
^  + 

ture  cf  scale  sizes  of  tens  of  km.  The  O  is  decreased  in  the  Fg  peak  due  to  the 
effect  of  the  electric  fields  that  increase  the  loss  rate  of  the  ion  molecule  reactions 
with  Nj  and  O2.  The  substantial  increase  in  the  N2^  is  due  to  enhanced  ionization 
by  the  precipitating  particles. 

3.4  Summary 

Measurements  of  the  total  density  below  90  km  reveal  substantial  structure 
with  scale  sizes  that  are  in  agreement  with  theoretical  predictions  of  gravity  wave 
theory  and  tidal  theory.  At  higher  altitudes  gravity  waves  are  also  apparent  but 
with  longer  wavelengths  (tens  of  km).  Satellite  observations  reveal  wave  structure 
in  the  species  concentrations  with  scale  sizes  from  tens  of  km  to  several  hundred 

45.  Trinks,  H. ,  Mayr,  H.G.  and  Philbrick,  C.R.  ( 197 8a >  Momentum  source 

signatures  in  thermospheric  neutral  composition,  J.  Geophys.  Res. 
^;1641-1646.  - 

46.  Trinks,  H.  ,  Offerman,  D. ,  von  Zahn,  U.  and  Steinhauer,  C.  (1978b)  Neutral 

composition  measurements  between  90  and  220  km  altitude  by  rocket-borne 
mass  spectrometer,  J.  Geophys.  Res.  83:2169-217  6. 
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Satellite  Flight 


kilometers.  The  description  of  these  small-scale  features  of  atmospheric  varia¬ 
tions  would  require  a  great  deal  of  knowledge  concerning  the  sources  of  the  parti¬ 
cular  disturbance,  as  well  as  detailed  3-D  calculations  of  atmospheric  motion. 

This  is  not  feasible  at  the  present  time  and  one  must  roly  on  a  statistical  descrip¬ 
tion  of  the  magnitude,  spatial  scale-size  distribution,  and  frequency.  W  ith  such 
a  description  it  should  be  possible  to  evaluate  the  variability  of  the  resultant  in¬ 
frared  radiance. 

Figures  17  and  18  give  indications  of  the  scales  and  amplitudes  of  waves  and 
other  structure  features  that  have  been  observed  in  the  upper  atmosphere.  In 
Figure  17  the  range  of  theoretically  predicted  scale  sizes  for  gravity  waves  (Hines, 
1960^^)  and  tidal  oscillations  (Forbes,  1978^)  are  shown  for  both  the  vertical  and 
horizontal  structure.  On  the  outside  edge  of  the  figure,  several  of  the  structure 
scales  that  have  been  measured  (most  are  from  examples  shown  earlier  in  this 
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interesting  point  to  notice  is  that  at  low  altituaes.  itie  eliai’ai'teristie  .-cale  sivi’s 
■Ti’e  typienilv  oik  to  a  li  v  km  and  amnlitufies  are  a  lew  pi  r.'em.  !■  i  '..ivine  lo 
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Figure  17.  Theoretical  Scale  Sizes  for  Oravitv 
Waves  and  Tidal  Oscillations.  Outside  engcs 
indicate  measured  scale  sizes. 
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higher  altitudes  the  amplitude  and  scale  size  both  increase.  While  it  is  certainly 
not  a  complete  picture,  it  does  present  an  idea  of  the  tyiiical  behavior  undet 
moderately  disturbed  atmospheric  conditions.  These  considerations  apply  to 
the  total  density  of  the  atmosphere.  The  structure  and  variability  of  an  individual 
minor  species  may  be  substantially  greater. 

From  this  review,  we  conclude  that  substantial  variations  occur  iti  the  den¬ 
sity  distribution  of  the  infrared  radiating  species  of  the  atmosphere,  both  spatially 
and  temporally.  This  chapter  has  only  discussed  the  variations  and  scales  asso¬ 
ciated  with  the  species  concentrations,  but  the  excitation  processes  may  be  more 
important  to  the  radiance  under  some  circumstances.  At  present  it  is  difficult 
to  determine  the  full  impact  of  these  variations  on  the  resultant  infrared  radiance. 
It  can  be  inferred,  for  example,  that  the  density  of  CO2  at  150  km  would  be  ex¬ 
pected  to  increase  by  more  than  a  factor  of  ten  under  high  latitude  geomagnetic 
heating.  With  the  calculations  discussed  elsewhere  in  this  report  and  with  the  data 
presented  here  these  questions  can  begin  to  be  answered. 
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Figure  18.  Altitude  Distribution  of  Measured  Scale  Sizes  from  Indicated  Sources. 
Amplitude  of  structures  are  indicated  by  length  of  plotted  line.  The  rocketborne 
data  show  vertical  structure  and  the  satellite  shows  horizontal  structure 


47.  Philbrick,  C.R.,  Faire,  A.E.,  and  Fryklund,  D.H.  (1978a)  Measurements 

of  Atmospheric  Density  at  Kwajalein  Atoll,  18  May  1977.  AFGL-TR-0058 
AD  Ao54  784.  - - - - 

48.  Philbrick,  C.  R.  (1976)  Recent  satellite  measurements  of  upper  atmospheric 

composition.  Space  Research  XVI:289-295. 
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4.  ATMOSPHERIC  BACKGROUND  CLUTTER  ANALYSIS 


50 

Preliminary  calculations  (Titus  and  Anapol,  1979  )  of  the  noise  or  clutter 

produced  by  fluctuations  in  the  CO2  radiation  have  been  made  using  a  model  based 

on  satellite  mass  spectrometer  measurements  of  density  variations  in  the  upper 

atmosphere  combined  with  estimates  of  the  mean  radiance  along  a  path  obtained 

from  the  High  Altitude  Radiance  Model.  CO^  is  the  dominan  adiator  throughout 

much  of  the  infrared,  but  it  is  one  of  the  most  difficult  molecules  to  measure 

using  mass  spectroscopy.  Because  of  this  difficulty,  CO2  dt-nf  ity  data  at  higher 

altitudes  are  almost  nonexistent  and,  therefore,  Ar  data  (Figure  19(a))  was  taken 

to  be  approximately  representative  of  CO„.  As  discussed  in  Section  3,  Reber  et  al 
43  ^ 

(1975)  have  published  spatial  density  variations  for  AR  at  high  altitudes  (see 
Figure  10).  Using  the  assumption  that  the  radiance  is  directly  proportional  to 
molecular  density,  the  power  spectrum  (PSD)  of  the  radiance  spatial  structure 
was  derived  from  these  data.  The  correlation  length  is  found  to  be  approximately 
10  km  and  the  roll-off  of  the  PSD  is  f  ^  (Figure  19(b))  where  f  is  the  spatial  fre¬ 
quency.  Satellite  ionosphere  density  (SRI  wideband  experiment)  measurements  at 
these  altitudes  and  accelerometer  and  calorimeter  measurements  appear  to  con¬ 
firm  these  values.  The  large-scale  correlation  length  of  the  AE-C  mass  spectro¬ 
meter  data  is  approximately  450  km.  This  length  is  corroborated  by  tlie  SPIRE 
data  and  S3-1  mass  spectrometer  measurements. 

The  impact  of  radiance  structure  on  a  sensor  inbedded  in  the  limb  is  as 
shown  in  Figure  19(c).  A  shell  model  of  radiance  is  constructed  at  correlation 
intervals  and  the  power  spectrum  calculated  for  each  intersection  of  the  ;ine-of- 
sight  and  a  shell.  The  power  spectrum  is  integrated  at  each  range  point  to  obtain 
the  mean  square  clutter  in  the  sensor  passband  and  these  values  summed  along 
range  to  obtain  the  integrated  mean  square  clutter  noise. 

Measurements  during  aurorally  active  conditions  exhibit  much  larger  density 
variability  and  smaller  correlation  lengths.  In  addition,  the  fluctuations  of  Ar 
were  observed  to  be  much  larger  than  those  of  N.  This  mass  dependence  trans¬ 
lates  into  a  Spectral  dependent  model  of  clutter  since  each  species  will  behave 
differently.  Large  radiance  enhancements  during  aurora  have  been  measured 
at  shorter  wavelengths  and  are  theoretically  expected  at  longer  wavelengths. 
Therefore  (and  not  unexpected),  the  impact  of  the  disturbed  upper  atmospheric 
conditions  will  be  more  severe  than  the  quiescent  atmosphere. 


50.  Titus,  J.  and  Anapol,  M,  (1979)  Private  communication. 


BACKGROUND  CLUTTER 


tmosphere  Background  PSD 

Figure  19.  Atmosphere  Background  Clutter  Analysis 


5.  AURORAL  EXaTATION  OF  INFRARED  RADIATORS 


Even  if  the  concentration  and  density  distribution  of  each  of  the  potential  in¬ 
frared  emitters  were  known,  the  radiance  at  high  latitudes  may  be  difficult  to  pre¬ 
dict  because  of  the  fluorescence  and  direct  collisional  excitation  sources  together 
with  localized  temperature  enhancements.  In  this  section  a  brief  discussion  will 
be  given  of  several  sets  of  information  that  are  related  to  the  energy  sources  that 
input  to  the  upper  atmosphere.  Tlie  energy  input  in  the  auroral  zone  may  be  studied 

by  examining  the  ionospheric  response,  particularly  the  N„^  enhancement  due  to 

+  ^ 

particle  ionization  and  the  O  depletions  due  to  the  electric  field  produced  increases 
in  ion  molecule  loss  reactions.  The  visible  radiation  from  all-sky  cameras  on 
DMSP  satellite  photographs  provide  an  excellent  means  of  characterizing  the  loca¬ 
tion,  structure,  and  intensity  of  energy  deposition.  Ground-based  video  recording 
systems  have  also  been  used  to  study  the  auroral  characteristics.  A  recent  satel- 

51 

lite  ultraviolet  photometer,  Huffman  et  al,  (1979),  has  also  provided  nadir  look¬ 
ing  data  on  the  ultraviolet  emission  excited  in  the  auroral  zones.  In  Figure  20 
measurements  of  the  radiation  at  155  nm  (the  Lyman-Birge-Hopfield  band  of  mol¬ 
ecular  nitrogen),  from  a  20°  portion  of  the  auroral  oval  showing  arcs  and  contin¬ 
uous  aurora,  are  presented. 

An  aurora  is  a  region  of  considerable  variability  of  luminosity.  In  spatial  ex¬ 
tent  it  may  vary  in  width  from  the  order  of  km  to  hundreds  of  km,  and  in  length  in 
thousands  of  km.  The  structure  is  very  heterogeneous  and,  thus,  the  aurora  is 
potentially  a  significant  source  of  fluctuating  infrared  radiation.  Examples  which 
show  the  variety  of  auroral  forms  that  may  be  encountered  are  shown  in  Figure  2  1, 
from  the  DMSP  satellite  system  operating  in  the  infrared.  The  land  mass  outlined 
by  the  lights  is  the  eastern  part  of  the  United  States. 

In  addition  to  the  well-known  visible  aurora,  typified  by  the  DMSP  photographs, 
we  may  have  the  occurrence  of  midlatitude  auroral  red  arcs  during  periods  of  high 

magnetic  activity.  These  occur  to  the  south  of  the  southern  extension  of  the  visible 
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aurora.  Figure  22,  (N.  B.S. ,  1962'  )  have  widths  of  the  order  of  several  hundred 
km,  and  extend  for  thousands  of  km  in  longitude.  These  arcs  are  characterized  by 
a  predominance  of  the  630  nm  O  line.  Their  altitudes  are  typically  above  300  km. 

In  Figure  23  ion  concentrations  are  shown  for  two  satellite  flights  through  auro¬ 
ral  disturbances.  In  both  cases  the  aurora  were  Class  III.  The  ion  concentrations 
display  the  high  degree  of  irregularity  to  be  expected  under  such  circumstances. 


51.  Huffman,  R.E.,  LeBlanc,  F.  J. ,  Larrabee,  J.  C.  and  Paulsen.  D.E.  (1979) 

Satellite  Atmospheric  Radiance  Measurements  in  the  Vacuum  Ultraviolet. 
AFCL-TR -79-0151,  ADA078  628. 

52.  N.B.S.  (1962)  Auroral  arcs  in  mid-latitudes.  National  Bureau  of  Standards 

Technical  News  Bulletin  46:42-46. 
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Figure  20.  Radiation  at  155  nm  (Lyman-Birge-Hopfield)  for  a  20°  Portion  of  the 
Auroral  Oval 

5.1  Application  of  VuiMe  Monochromatic  Auroral  Imaging  Data  for  Modeling  infrared  Radiance 

Effective  design  of  sensitive  infrared  sensors  operating  in  the  limb  viewing 
mode  requires  a  statistically  significant  data  base  on  the  auroral  infrared  earth- 
limb  background.  Research  that  has  been  performed  under  the  iNA/AFGL  auroral 
measurements  program  has  shown  that  there  are  strong  auroral  emissions  near 
2. 7  and  4.  3  nm  (Stair  et  al,  1981^).  It  is  possible  to  develop,  to  first  order,  a 
comprehensive  2. 7  and  4.  3  nm  earthlimb  model  that  is  based  on  an  existing  all¬ 
sky  TV  data  base.  The  TV  data  base  includes  approximately  200  hours  of  all-sky 
multispectral  image  intensified  TV  imaging  data  of  aurora  and  airglow  that  have 
been  taken  and  stored  on  video  tape  in  the  time  period  1972  through  1977.  Methods 
for  modeling  (blue  redjnput  model,  BRIM)  infrared  emissions  at  2.  7  and  4.  3  nm 
on  the  basis  of  ground-based  photometric  measurements  of  the  aurora  in  the  spec¬ 
tral  regions  blue  (4278A)  and  red  (6300A)  have  recently  been  developed  and  verified 
against  Infrared  data  that  were  obtained  simultaneously  via  rocketborne  sensors. 
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Figure  21.  DMSP Infrared  Photographs  ofVarious  Auroral  Forms 


Figure  22,  The  Red  Arc  (shaded  region) 

Lies  South  of  the  Southern  Extension  of  the 
Visible  Aurora  (_) 

In  this  section  the  use  of  the  BRIM  capability  to  construct  images  of  the  infra¬ 
red  earthlimb  radiance  from  successive  pairs  of  blue  and  red  all-sky  TV  images 
is  discussed.  The  images  are  obtained  at  a  rate  of  once  every  15  sec  or  faster. 
The  spatial  resolution  of  the  TV  is  such  that  structure  on  the  scale  of  several  km 
may  be  resolved  in  the  radiance  images.  The  all-sky  TV  coverage  has  a  zenith 
full  cone  angle  of  160°.  This  is  extensive  enough  so  that  images  of  horizontal  ex¬ 
tent  of  up  to  7  00  km  and  tangent  altitude  from  40  to  140  km  may  be  constructed. 
Spatial  power  spectral  densities  (PSDs)  may  be  obtained  from  the  images  and  tem¬ 
poral  PSDs  may  be  obtained  from  a  series  of  images. 

53 

Methods  inferred  from  the  work  reported  by  Rees  and  Luckey  (1974)  may  be 
used  to  approximately  construct  an  energy  deposition  profile  q(z,t)  from  the  TV 
brightness  data  measured  in  narrow  spectral  bands  in  the  blue  (4278A)  and 
in  the  red  (6300A)  spectral  region.  The  "blue"  data  4jrljj(t)  are  proportional  to  the 
height  integrated  energy  deposition  and  the  ratio  Ig/Ip  of  the  blue  to  red  may  be 
associated  with  the  altitude  where  q  is  a  maximum  if  one  assumes  that  the  energy 
dependence  of  the  incident  electrons  takes  the  form  <t>  =  '^’pEe  ''  .  Kumer  (1977)' 

53.  Rees,  M.H,  and  Luckey,  D.  (1974)  Auroral  electron  energy  derived  from 

ratio  of  spectroscopic  emissions,  J.  Geophys.  Res.  79:5181-5186. 

54.  Kumer,  J.  B.  (1977)  Further  Eva^ation  of  ICECAP  Auroral  4.  3  pm  Zenith 

Radiance,  Haes  Report  No.  57  DNA  4260F. 
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discussed  in  detail  the  BRIM  method  for  approximately  constructing  a  deposition 
profile  q(z,t)  from  the  blue  and  red  brightness  data.  Examples  of  the  excellent 
agreements  between  the  up  and  down  leg  BRIM  calculations  for  4.  3  txm  zenith  radi¬ 
ance  and  the  data  obtained  in  the  case  of  the  an  auroral  event  of  25  February  1974 
are  shown  on  Figures  24  and  2  5.  Ground-based  photometric  red  and  blue  auroral 
data  taken  at  the  point  of  up-leg  and  downleg  penetration  at  80  to  110  km  are  input 
to  the  BRIM  model. 

The  goal  is  to  use  a  series  of  TV  frames  of  blue  data  4irlg  to  generate  a  dyna¬ 
mic  image  of  the  corresponding  infrared  earthlimb  radiance  as  it  would  appear  when 
viewed  from  a  satellite  at  some  considerable  distance  from  the  earth.  Line-by-line 
(or  selected  line)  digitization  of  auroral  images  and  image  sequences  is  accomplished 
with  a  system  which  incorporates  a  controller  that  selects  the  proper  TV  frame,  a 
digitizer,  and  an  output  couple  to  transfer  the  digital  data  stream  to  digital  magnetic 
tape. 


Figure  24.  25  February  1974  Upleg  4.3>im  Zen¬ 
ith  Radiance  Data  (curve  D)  are  Compared  With 
Curves  Ay,  the  BRIM  Calculation  for  the  Auro¬ 
ral  Component  of  Zenith  Radiance;  OH.  the  Cal¬ 
culated  Component  Due  to  Vibration  Transfer 
from  OH"*";  Ayj(xl0),  the  Ambient  Thermal  Com¬ 
ponent  Multiplied  by  a  Factor  10  So  That  It  is  on 
Scale;  and  T^,  the  Sum  OH  +  A^  +  A^^ 
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Figure  25.  A  Comparison  of  the  2  5 
February  1974  Downleg  Data  (curve  D) 
and  Calculations  (curve  BRIM)  Based 
on  Our  BRIM  Model 

To  date  Kumer  and  associates  have  developed  software  for  the  calculation  of 
2.7  and  4.  3  Pm  earthlimb  radiance  from  digitized  blue  and  red  zenith  viewing  all¬ 
sky  TV  data.  Due  to  the  lack  of  mass  digitization  of  the  TV  data  base  at  this  time 
a  version  of  the  simulated  TV  data  has  been  utilized  to  test  the  software.  In  this 
simulation  of  the  earthlimb,  radiance  is  calculated  along  51  lines-of-sight  with 
tangent  altitude  110  km.  The  middle  line-of-sight  passes  directly  over  the  zenith 
of  the  TV  station.  The  surface  projection  of  this  line-of-sight  runs  from  magnetic 
northwest  to  southeast  and  is  30°  magnetic  south  of  due  magnetic  east.  The  lines- 
of-sight  are  parallel,  2  km  apart,  and  tangent  at  110  km  altitude. 

The  earthlimb  radiance  calculations  are  performed  on  50  pairs  of  simulated 
blue/red  digitized  TV  data  frames.  These  pairs  of  frames  are  spaced  in  time  by 
2  0  sec  and  thereby  define  an  interval  0  <  t  <  16.  67  min.  The  first  pair  was  gener¬ 
ated  from  scanning  photometer  data.  This  initial  pair  is  considered  to  represent 
a  steady-state  condition  that  had  persisted  for  many  min  prior  to  t  =  0  min. 

At  t  =  0  it  is  assumed  that  a  multiple  arc  structure  appears  in  the  field  of  view 
(FOV)  and  commences  to  move  across  the  FOV  in  the  positive  direction  perpendicular 
to  the  lines-of-sight.  In  the  blue  and  red  TV  data  the  multiple  arc  structure  is  re¬ 
presented  sinusoidally  by 


10  (1  -  cos  2a’k(x  '  vt}) 


4»Ib  = 


and 


4irlB  =  4)rlg/l0  , 

where  k  =  0.  1  km  =  7.  14  km/min,  and  the  brightness  units  are  kilo -Rayleighs 
(kR).  The  structure  is  assumed  uniform  along  lines  parallel  to  the  lines-of-sight. 

Results  for  the  earthlimb  radiances  calculated  on  the  middle  line-of-sight 
(that  is,  X  =  0)  on  the  time  interval  0  <  t  <  16.67  min  are  shown  on  Figure  26.  One 
can  see  the  2.7  pm  emission  tracks  the  temporal  structure  of  the  simulated  moving 
auroral  arcs.  This  results  from  the  fact  that  the  response  time  of  the  N(  D)  mech¬ 
anism  for  2.  7  pm  auroral  emission,  about  0.  5  sec  at  210  km,  is  considerably 
shorter  than  the  period  (-1.  67  min)  of  the  simulated  aurora.  The  response  time 
for  the  C'O.,  4.3  pm  aurora  at  110  km  is  approximately  16  min,  considerably  longer 
than  the  period  of  the  auroral  excitation.  This  explains  the  temporal  behavior  of 
the  auroral  CO„  4.  3  pm  earthlimb  radiance  as  shown  on  Figure  26.  This  radiance 
asymptotically  approaches  an  equilibrium  DC  level  (si.  8  X  10  w/ cm“  sr)  near 
t  =  16.  67  min.  It  also  displays  a  small  1.  67  min  period  ripple  in  response  to  the 
aurora.  The  amplitude  of  the  COj  4.  3  pm  ripple  is  about  30  times  less  than  the 
amplitude  of  the  NO  2.  7  pm  radiance.  Also  despite  the  fact  that  it  is  well  known 
that  the  DC  level  for  vertically  viewed  auroral  COg  4,  3  pm  emission  is  much  lar¬ 
ger  than  for  NO  auroral  2.  7  pm,  we  see  that  in  viewing  on  the  limb  CO2  4.  3  pm 
self  absorption  renders  the  average  radiance  values  about  equal  for  these  two 
emissions. 

The  steady-state  spatial  structure  at  t  =  o  is  shown  on  Figure  27.  The  posi¬ 
tive  gradient  results  from  the  north  sky  position  of  measurements  of  an  arc  of 
26  October  1978  with  respect  to  the  position  of  the  scanning  photometer  that  ob¬ 
tained  the  data  that  are  used  to  generate  the  first  frame  pair  of  simulated  digitized 
TV  data. 
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TEMPORAL  STRUCTURE 


min 


Figure  26.  Temporal  Dependence  on  Line -of-Siglit  Passing 
Directly  Overhead  TV  Station  at  110  km  Tangent  .Mtitude  on 
the  0  £  t  5  16.  63  min  Interval.  Projected  view  direction  i.® 
approximately  12 oo  southeast  of  due  magnetic  north 
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Figure  27.  Horizontal  Structure  in  Simulated  Steady-State  at 
t  =  0,  Based  on  Photometric  Data  Obtained  26  October  1978 
by  R.  D.  Sears^^ 
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55.  Scars,  R.D.  (1979)  Private  communication. 
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6.  Summary  and  Recommendations 


In  this  review  of  atmospheric  variability  and  its  effect  on  the  variability  of 
atmospheric  infrared  radiation,  certain  important  problems  are  evident.  The  i 

most  important  problem  is  the  basic  lack  of  measurements  concerning  the  vari¬ 
ability  of  the  infrared  radiance  from  the  upper  atmosphere.  Such  measurements 
are  highly  desirable  and  should  be  made  at  the  earliest  opportunity.  Until  these 
measurements  are  made  the  prediction  of  the  atmospheric  infrared  radiation  must 
rely  on  numerical  codes  such  as  the  High  Altitude  Infrared  Radiance  Model  de¬ 
scribed  in  Chapter  2.  This  model  is  well  established  and  predictions  made  using 
it  are  in  general  agreement  with  the  limited  measurements  that  have  been  per¬ 
formed.  However,  the  calculations  do  not  include  the  effects  of  atmospheric 
variability.  This  must  be  included  in  future  predictions  as  the  capability  of  most 
systems  of  interest  is  limited  principally  by  variations  in  the  radiance  that  appear 
as  false  signals. 

The  models  that  describe  the  long-term  variations  in  the  atmospheric  proper¬ 
ties,  which  are  the  necessary  input  to  codes  that  predict  atmospheric  radiation, 
are  discussed  in  Chapter  3.  The  predictions  of  these  models  are  in  satisfactory 
agreement  with  measurements  for  the  long-term  variations;  however,  the  models 
do  not  attempt  to  follow  the  short-term  variations  such  as  would  be  seen  in  frame - 
to-frame  viewing.  In  particular,  the  Cole-Kantor  low  altitude  model  and  the  J77 
model  seem  quite  satisfactory  for  the  prediction  of  the  long-term  variability. 

However,  a  problem  of  how  to  Join  these  models  near  the  homopause  remains  and 
must  be  examined  further  before  satisfactory  predictions  of  the  radiance  can  be 
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made.  The  joining  of  the  properties  of  these  two  models  is  currently  under  study. 
As  both  of  these  models  describe  only  average  or  mean  conditions  ovei'  long  peri¬ 
ods  of  time,  the  most  serious  problems  that  arise  in  the  use  of  disjoint  models  is 
in  the  calculation  of  spatially  integrated  radiance  values  and  the  interpretation  of 
the  calculated  results  of  radiance  when  the  limb  passes  through  that  transition 
region. 

A  second  problem  in  the  application  of  the  models  arises  in  prediction  of  the 
long-term  variability  of  the  minor  species  that  are  the  major  contributors  to  the 
infrared  radiation.  These  predictions  are  based  principally  on  theoretical  models 
of  atmospheric  behavior  and  a  limited  number  oi  measurements.  However,  these 
models  are  questionable  for  the  minor  species  and  must  be  examined  in  greater 
detail.  Preliminary  predictions  of  long  term  variability  of  the  infrared  radiation 
that  are  presented  in  Chapter  4  show  ihe  degree  to  which  the  limb  radiance  varies 
with  altitude.  The  application  of  these  (  Its  would  change  the  effective  limb 
altitude  of  a  system  that  makes  mcssureni.cnts  near  the  horizon. 

Data  do  e-;ist  on  the  long-term  variability  of  minor  species  but  at  present  they 
have  not  been  fully  analyzed.  The  qu.<».ntity  of  data  is  still  not  sufficient  to  allow 
the  minor  species  to  be  modeleu  to  the  degree  that  the  major  components  of  the 
atmosphere  are  modeled.  At  the  lower  altitudes  these  densities  are  best  deter¬ 
mined  from  rocket  flights  that  should  be  continued  in  order  to  collect  the  amount 
of  data  of  the  time  and  space  that  is  required.  At  higher  altitudes  satellite  data 
are  available  and  are  being  analyzed.  However,  most  of  the  satellite  data  was 
obtained  at  altitudes  above  those  where  the  principle  radiation  is  emitted;  thus,  the 
total  amount  of  data  is  limited.  Under  certain  conditions,  means  other  than  direct 
sampling,  such  as  radar  returns  and  optical  measurements,  are  available  and 
should  be  pursued.  Theoretical  models  that  describe  the  chemistry  and  transport 
are  also  available  and  can  be  used  under  certain  circumstances. 

Another  serious  problem  is  the  degree  of  the  short-term  temporal  variability 
that  would  affect  the  frame-to-frame  operation  of  staring  measurement  systems 
and  the  small-scale  spatial  variability  that  would  affect  the  operation  of  scanning 
measurements.  The  available  data  concerning  species  concentrations  are  dis¬ 
cussed  in  Chapter  5  and  preliminary  conclusions  as  the  expected  degree  of  the 
variability  presented.  There  are  data,  such  as  the  AU-C,  -D,  -E,  and  S3- 1 
satellite  measurements,  that  should  be  analyzed  further.  The  preliminary  results 
of  this  analysis  were  used  to  construct  a  simplified  model  of  the  COg  radiance  as 
discussed  in  Chapter  5  and  applied  to  the  calculation  of  a  spatially  integrated  power 
spectral  density  for  analyses  of  certain  systems.  It  must  be  emphasized  that  the 
data  presently  available  concerning  the  short-term  variability  are  those  on  the 
species  concentration  and  not  on  the  infrared  radiance  variability.  Thus,  data  on 
the  short-term  variability  of  the  infrared  radiance  is  very  important  before  the 
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predictions  of  the  radiance  from  species  concentration  can  be  fully  accepted. 

These  predictions  do  rely  on  basic  knowledge  of  the  processes  that  lead  to  infrared 
radiation.  This  study  must  be  pursued  and  the  results  folded  into  the  predictive 
models.  Of  course,  this  process  requires  measurements  to  be  made  of  both  the 
source  and  the  emitted  radiation,  but  not  simultaneously  for  all  cases. 

Again,  in  the  determination  of  the  short-term  variability  many  types  of  meas¬ 
urements  can  be  used.  The  most  useful  source  of  data  on  the  atmospheric  con¬ 
stituents  are  measurements  of  species  concentrations  of  the  infrared  emitters  on 
a  spatial  and  temporal  scale  that  overlaps  with  system  measurement  parameters. 
Spatial  dimensions  are  presently  adequate  and.  as  more  data  are  available  from 
satellite  and  rocket  flights,  satisfactory  models  should  be  forthcoming.  Temporal 
measurements  for  time-scales  of  interest  must  be  made  from  fixed  platforms 
(such  as  the  earth,  balloons,  and  possibly,  the  shuttle)  using  inferential  data  such 
as  radar  returns  or  optical  measurements.  Such  optical  measurements  can  be 
made  at  wavelengths  other  than  those  of  current  system  interest.  Incoherent 
radar  returns  can  also  yield  the  atmospheric  temperature  that  is  a  crucial  para¬ 
meter  in  the  prediction  of  the  infrared  radiance. 

For  naturally  disturbed  atmospheric  conditions,  such  as  an  aurora,  coordi¬ 
nated  measurements  are  available  and  more  are  being  planned.  These  measure¬ 
ments  involve  both  data  on  the  energy  input  and  the  resulting  radiation,  as  well  as 
measurements  that  can  be  used  to  infer  local  conditions,  such  as  electron  densities, 
within  the  radiation  source  region.  In  Chapter  5  a  method  of  using  some  of  these 
data  to  predict  the  spatial  and  temporal  dependence  of  the  limb  radiance  at  2.  7  and 
4.  3  Mm  is  discussed. 

This  summary  and  recommendations  of  possible  measurements  must  be  recog¬ 
nized  as  being  based  on  incomplete  information.  As  new  measurements  are  made 
and  their  results  added  to  our  knowledge  of  atmospheric  radiances,  new  routes  of 
investigation  may  be  discovered. 
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